Update: Precision D s decay constant from full lattice QCD using very fine lattices 
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We update our previous determination of both the decay constant and the mass of the D s meson 
using the Highly Improved Staggered Quark formalism. We include additional results at two finer 
values of the lattice spacing along with improved determinations of the lattice spacing and improved 
tuning of the charm and strange quark masses. We obtain mo s = 1.9691(32) GeV, in good agreement 
with experiment, and fo e = 0.2480(25) GeV. Our result for fo B is 1.6a lower than the most recent 
experimental average determined from the D s leptonic decay rate and using V C3 from CKM unitarity. 
Combining our fr> B with the experimental rate we obtain a direct determination of V C3 = 1.010(22), 
or alternatively 0.990 ^q qib using a probability distribution for statistical errors for this quantity 
which vanishes above 1. We also include an accurate prediction of the decay constant of the n c , f Vc 
= 0.3947(24) GeV, as a calibration point for other lattice calculations. 



I. INTRODUCTION 

Lattice QCD is now a firmly established method for 
providing precision tests of the Standard Model [T] . Com- 
bined with experiment, lattice QCD calculations have the 
potential to uncover new physics provided that both the 
theoretical and experimental results are accurate enough. 

The most accurate lattice QCD calculations are those 
for the masses of 'gold-plated' mesons, where few MeV er- 
rors are now possible across the entire spectrum. This ac- 
curacy is at the level where electromagnetic effects on the 
meson masses, currently missing from lattice QCD calcu- 
lations, have to be estimated and included. Ref. [2] gives 
a recent summary including predictions of masses that 
have been made ahead of experiment. The meson masses 
are extracted from simple 'two-point' hadron correlation 
functions calculated on the lattice from combining ap- 
propriate valence quark and antiquark propagators. An- 
other parallel set of quantities that can be determined 
from the same correlation functions are the meson de- 
cay constants. Calculations of these can be compared to 
experimental results for rates of annihilation to photons 
for neutral unflavored vector mesons and to W bosons for 
charged pseudoscalars. By determining as complete and 
accurate a picture as possible for decay constants along 
with masses we provide a stringent test of the Standard 
Model. Physics beyond the Standard Model can intro- 
duce new ways to decay to leptons for some mesons, and 
so accurate comparison of decay constants between the- 
ory and experiment can also provide direct constraints 
on new physics models. 
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Here we focus on results for one quantity, the decay 
constant of the D s meson, fo B , which has been a show- 
case for the impact that accurate lattice QCD calcula- 
tions can have, particularly when ahead of experimental 
results. We will update our result from 2007 [3], mak- 
ing several improvements to the calculation. It is impor- 
tant to understand that /rj s is not calculated in isolation; 
as discussed above, it is one piece of the range of QCD 
physics that is calculable on the lattice. The other pieces, 
where they can also be tested against experiment, lend 
weight to the confidence we have in our error analysis. 
This is particularly true for our calculation because we 
can calculate a range of different quantities all with the 
same method. So here we also update our results for the 
mass of the D s meson and discuss other calculations that 
will provide further tests. First we review briefly some 
background to the calculation of fn s . 

Decay constants for light pseudoscalar mesons (/ w and 
/k) have been calculable with errors at the few percent 
level since 2004 [3J, being one of the first calculations 
done in lattice QCD once ensembles of gluon field config- 
urations were available that included the full effect of u, 
d and s sea quarks with a light enough mass for the u/d 
quark to enable controlled extrapolation to the physical 
point. These calculations were done using the improved 
staggered (asqtad) formalism [5J [B] which has a number 
of advantages over previous formalisms, that mean that 
the calculation of f„ and fx can be done accurately. Key 
requirements for these calculations are a quark formalism 
(such as improved staggered quarks) which : 

• has an absolutely normalised operator to couple to 
the W boson; 

• is improved so that it has small discretisation errors 
(0(cv s a 2 ) for improved staggered quarks) and 

• is numerically fast so that large ensembles of gluon 
field configurations can be made including sea 
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quarks and so that many hadron correlation func- 
tions can be calculated per configuration, for small 
statistical errors. In addition a large physical vol- 
ume (> 2.5fm 3 ) is necessary so that finite volume 
effects arc reduced to the 1% level. Having all of 
these features means that accurate extrapolations 
to the physical continuum limit can be made. 

Calculations of decay constants for mesons containing 
the heavier c quark became important with the promise 
of results from the CLEO-c experiment. The first lattice 
results for fjj and fn s appeared from the Fermilab Lat- 
tice/MILC collaborations in 2005 as predictions ahead of 
experiment [71. They used the 'Fermilab' formalism [8], 
developed many years previously for heavy quark physics, 
and had errors of 8%. This led to the unfortunate im- 
pression that decay constants for D and D s mesons were 
inevitably much less accurate than those for ir and K and 
errors would only be slowly reduced as higher statistics 
and the advent of finer lattices reduced statistical errors 
and systematic errors from discretisation effects. Because 
the Fermilab formalism predated the improved staggered 
formalism, however, these calculations had not made use 
of any of the features discussed above that made f v and 
Jk so accurate. 

For c quarks the issue of discretisation errors becomes 
more important than for the lighter quarks. In 2007 we 
showed that further improving the improved staggered 
formalism to the Highly Improved Staggered (HISQ) for- 
malism [3] produces a quark formalism that has all the 
good features of the asqtad formalism outlined above but 
also significantly smaller discretisation errors. In fact the 
discretisation errors are small enough that HISQ can be 
used for c quarks as well as u/d and s quarks and using 
the same formalism for all 4 lightest quarks has enormous 
advantages. We used HISQ for all the valence quarks to 
calculate all 4 decay constants: f v , f K , f D and f Ds to 
better than 2% accuracy [5]. Our results were: 

U = 132(2)MeV (1) 

f K = 157(2)MeV 

f D = 207(4)MeV 

f Da = 241(3)MeV 

Although Jd and fr>, still have noticeably larger discreti- 
sation errors (and therefore contributions to the system- 
atic error from the extrapolation to the a — > limit) than 
fx and /,r there are smaller systematic errors from, for 
example, finite volume effects. This leads then to the 
expectation, and the result, of very similar final errors. 
Our error for fjj a was somewhat smaller than that for 
Jd (1.3% versus 1.8%) since the D s contains no valence 
u/d quarks and is therefore much less sensitive to the chi- 
ral extrapolation to the physical u/d quark mass. This 
makes fo s a particularly accurate quantity to calculate 
in lattice QCD. 

Since, at that time, fjj and fjj s were only known to 
6-8% from experiment [TUHT2], we had the added test, 



unavailable to the Fermilab formalism, of agreeing with 
experiment for f n and Jk- An additional very stringent 
test that had not previously been done was the determi- 
nation of the mass of the D s and D mesons along with 
their decay constants. The masses are known to better 
than 1 MeV experimentally. We were able to achieve er- 
rors from lattice QCD of 7 MeV (0.3%) by determining 
the difference between the D or D s mass and one half 
that of the rj c . Electromagnetic effects on the masses, 
missing from the lattice QCD calculation, had to be al- 
lowed for in achieving this accuracy. Good agreement 
between lattice QCD and the experimental results was 
obtained. We quoted mn s — 1.962(6) GeV and vnu = 
1.868(7) GeV 0. 

Following our result much improved experimental 
numbers for f D [T3] (206(9) MeV) and f Ds (274(11) 
MeV) became available from CLEO |14) . This produced 
the exciting picture in the summer of 2008 that agree- 
ment between experiment and our result for fjj was very 
good but that the experimental result for fo s (includ- 
ing averages with results from BaBar and Belle [15]) was 
significantly larger than our lattice QCD value, see, for 
example, |16l 117] . Since the experimental errors were 
still much bigger than ours the discrepancy, of 3er, was 
dominated by the experimental error. A burst of activ- 
ity from other lattice QCD calculations produced results 
that agreed with ours but, having errors several times 
larger, often also agreed with the experimental one [IS] . 
This led to much speculation about the existence of new 
physics (that had to affect D s but not D) as well as 
limits on new physics from the fact that the experimen- 
tal fo e was larger (and not smaller) than the Standard 
Model result from our calculation [20] . 

Since then improved statistics and further results from 
other channels [2TH21] have brought down the experi- 
mental average and reduced its error to 2%. In early 
2010, the Heavy Flavor Averaging Group (HFAG) gave a 
world average result from experiment of fn s = 0.2546(59) 
GeV [25], 2cr above our 2007 result. In the meantime we 
have extended our lattice QCD calculation using HISQ 
quarks to even finer lattices as well as improving the accu- 
racy with which we determine the lattice spacing (which 
provides the calibration of the energy scale) and fix the 
c and s quark masses. This has led our result for fu a 
to move upwards, as we show here, to 0.2480(25) GeV, 
with a slight improvement in the error to 1%. The main 
reason for the upward shift is the recalibration of the 
lattice spacing. The experimental average as of Octo- 
ber 2010 has moved up again slightly with new results 
from BaBar [25]. Our value for fu s is now 1.6er from the 
experimental average and this reduces considerably the 
room for new physics in this quantity. 

In section [TT] we describe the lattice QCD calculation 
and in section Mil the results. These sections contain 
technical details which may not be of interest to those 
without a lattice QCD background. As well as fr>, we 
give results for mu a which, as discussed above, is an im- 
portant check on the calculation. We also show results 
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for f Vc , the decay constant of the rj c . This cannot be 
accessed directly from experiment but provides an ex- 
cellent 'figure of merit' for lattice QCD calculations in 
charm physics. We give the result to 0.6% so that other 
lattice QCD calculations can compare to this when quot- 
ing numbers for fjj a . In section IV we discuss the picture 



that emerges from the current experimental and lattice 
QCD results for fjj s , including the update we give here. 
We have tried to make this section readable by those 
that skipped the earlier technical details. We will also 
comment on the effects of the recalibration of the lattice 
energy scale on the other calculations included in [3J, i.e. 
Ik, fir an d fn- Section |V| gives our conclusions. 



II. LATTICE QCD CALCULATION 

We work with 11 different ensembles of gluon field con- 
figurations provided by the MILC collaboration. These 
include the effect of u, d and s sea quarks using the 
improved staggered (asqtad) formalism and the fourth 
root 'trick'. This procedure has passed various tests in- 
dicating that it is a valid discretization of QCD [28l430j . 
Configurations are available with large spatial volumes 
(> 2.4fm 3 ) for a wide range of values of the lattice spac- 
ing, a, and at multiple values of the sea light and strange 
quark masses. The u and d quark masses are taken to be 
equal in the sea (m u — m c \ = m{) for numerical speed. 
This has negligible effect on the calculations described 
here. We use configurations at 5 values of the lattice 
spacing between 0.15 fm and 0.045 fm with parameters 
as listed in Table H] We have chosen the ensembles so that 
we can test the dependence of our results on each of: the 
lattice spacing; the physical volume; the sea light quark 
mass and the sea strange quark mass. 

On these configurations we have calculated quark prop- 
agators for charm quarks and strange quarks using the 
HISQ action. The numerical speed of FIISQ means that 
we have been able to use several nearby quark masses 
for charm and strange to allow accurate interpolation to 
the correct values. This is described in the next section. 
These propagators are combined together to make pseu- 
doscalar meson correlators with valence quark content ei- 
ther 'charm-charm', 'charm-strange' or 'strange-strange'. 
By fitting the correlators as a function of the time sepa- 
ration of the source and the sink on the lattice we are 
able to determine the ground-state pseudoscalar mass 
(i.e. that of the ry c , D s or r) s ) and the amplitude with 
which the ground-state meson is created or destroyed by 
the local temporal axial current. This latter quantity is 
directly related to the decay constant. 

The HISQ action [9 is an extension of the asqtad im- 
proved staggered quark action, which is itself based on 
the unimproved (naive) staggered quark action. The 
unimproved staggered action is equivalent to a simple 
'naive' discretization of the continuum quark action to 



TABLE I: Ensembles (sets) of MILC configurations with size 
i 3 x T and sea mass parameters mf sq and mf sq used for this 
analysis. The sea ASQTAD quark masses (I = u/d) are given 
in the MILC convention where no is the plaquette tadpole 
parameter. Values of Mo are given in Table |VI| The lattice 
spacing values in units of ri after 'smoothing' are given in 
the second column [27]. Sets 1 and 2 are 'very coarse'; sets 
3, 4, 5, 6 and 7 'coarse'; sets 8 and 9 'fine'; set 10 'superfine' 
and set 11 'ultrafine'. The final column gives the number of 
configurations and the number of time sources per configura- 
tion used for calculating quark propagators for the best-tuned 
parameter sets on each ensemble. 



Set 


ri/a 


asq 

auom l 


au mf q 


L/T 


N cf x N t 


1 


2.152(5) 


0.0097 


0.0484 


16/48 


631 x 2 


2 


2.138(4) 


0.0194 


0.0484 


16/48 


631 x 2 


3 


2.647(3) 


0.005 


0.05 


24/64 


678 x 2 


4 


2.618(3) 


0.01 


0.05 


20/64 


595 x 2 


5 


2.618(3) 


0.01 


0.05 


28/64 


269 x 4 


6 


2.644(3) 


0.02 


0.05 


20/64 


600 x 2 


7 


2.658(3) 


0.01 


0.03 


20/64 


328 x 2 


8 


3.699(3) 


0.0062 


0.031 


28/96 


566 x 4 


9 


3.712(4) 


0.0124 


0.031 


28/96 


600 x 4 


10 


5.296(7) 


0.0036 


0.018 


48/144 


201 x 2 


11 


7.115(20) 


0.0028 


0.014 


64/192 


208 x 1 



give, on the lattice: 

S = (7 -A(17) + ma) ip(x). (2) 

X 

where ma is the quark mass in lattice units. A(U) is a 
discrete version of the covariant derivative coupling to the 
lattice gluon field U^x), which is a set of SU(3) matrices 
sitting on the links of the lattice: 



U l _ l {x)ijj(x + fx) 



fx)ip(x - /i)] ■ 
(3) 

In the improved staggered formalism the gluon field in 
the covariant derivative is smeared i.e. U^(x) is replaced 
by a sum of products of matrices tracing out more 
complicated paths between x and x + (i [5] . The smear- 
ing introduces a form factor that reduces the coupling 
between the quark and high momentum ( p ~ n / a ) glu- 
ons that cause a particular type of discretisation error 
for staggered quarks. This error in principle appears at 
a s a 2 but in practice is very large for unimproved stag- 
gered quarks. The error is seen most clearly in the mass 
differences between different 'tastes' of pseudoscalar me- 
son, created by different point-split pseudoscalar oper- 
ators. These mass splittings are proportional to a 2 and 
are strongly reduced on going from unimproved staggered 
quarks to improved staggered quarks [6j. Most smear- 
ing methods introduce additional discretisation errors. 
This is avoided here by the specific form of the smearing 
used 5 . In the highly improved staggered quark action 
this smearing is applied twice with a reunitarisation of 
the gluon field in between. We also apply a projection 
back on to SU(3) for the gluon field, although this makes 
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little difference in practice. We then find another further 
large reduction in the splittings between different tastes 
of pseudoscalar mesons [9]. In the pseudoscalar case the 
splitting in the squared masses (AmJ) is roughly con- 
stant (for quark masses that are not too large) and so the 
splittings in the pseudoscalar masses themselves (Am ff ) 
fall with quark mass. Thus these 'taste-changing' er- 
rors are generally smaller for charm quarks than strange 
quarks 131) , and they are particularly small with the 
HISQ action. 

Other, more mundane, discretisation errors are tackled 
by standard improvement techniques. A simple analy- 
sis in Fourier space of the symmetric difference of equa- 
tion [3] shows that this has errors of 0(a 2 ) which can be 
corrected by the addition of a (pa) 3 term. This term, 
known as the Naik term [32] . appears in the improved 
staggered quark action as a mixture of 3-link and 1-link 
differences. The improved staggered quark action then 
has discretisation errors that are 0(a 4 ), apart from ra- 
diatively generated errors at 0(a s a 2 ). The HISQ action 
uses the same Naik term (except that it contains smeared 
gluon fields) but corrects it further for discretisation er- 
rors when using quark masses appropriate to charm or 
heavier. Discretisation errors controlled by the quark 
mass ma become important in that case, and we adjust 
the coefficient of the Naik term so that it takes value 
(1 + e) instead of 1 [5]. Then, schematically 



S = ip(x) (7 • A(Z7) + ma) ip(x) 



where 



A„ 
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(4) 



(5) 



£ is a function of ma (starting at (ma) 2 ) calculated to 
give the correct quark dispersion relation (energy as a 
function of momentum) at tree level. Here we give an 
exact formula for e at tree level, etree, given an expansion 
for the tree level pole mass, m tree , as a function of the 
mass ma in the lattice action [5] : 



ma[l — 

1783 
537600 

4- 



80 (ma) + 2240 (ma) 



(6) 



(ma) s 



76943 
23654400 



(ma) w + . 



12mt ree a 



cosh(mt ree a) sinh(mt ree a) 

(sinh(m tree a)) 2 



(7) 



These equations are obtained by solving the condition 
for the 'kinetic mass', M2 = [d 2 E / dp 2 ]^ 1 , to be equal to 
the tree level pole mass, m tree . mtree in turn solves the 
pole condition at zero momentum. Including a Naik coef- 
ficient of (l+ttree) means that the leading (in the velocity 
expansion) (ma) 4 errors are removed in the HISQ case, 
and so remaining discretisation errors are suppressed ei- 
ther by a s or by the fact that heavy quarks are nonrela- 
tivistic in their bound states, e can be fixed nonpertur- 
batively by demanding that the 'speed of light' be 1, and 



this was done in earlier calculations [5]. However it was 
found that nonperturbative results for e were close to the 
tree level result in the HISQ case and so here we simply 
define e to take the value ttree above. 

It is numerically very fast to calculate quark propa- 
gators for staggered actions because they have only one 
spin component. This means that we can readily calcu- 
late propagators from several different time sources on 
the lattice for improved statistics. Table [T] lists the num- 
ber of configurations used from each ensemble and the 
number of time sources per configuration. To increase 
statistics further we use a 'random wall' source for the 
quark propagator instead of a delta function [4 . The 
random wall is a set of U(l) random numbers with unit 
norm on every point of the source time slice (separately 
for each color) and is used as the source for the inversion 
to calculate the quark propagator. The same random 
wall is used for all propagators from a given time source 
on a given configuration so that when any propagator is 
combined with the complex conjugate of another to form 
a meson correlator the random numbers cancel except 
where the initial spatial points and colors are the same. 
This effectively increases the number of meson correla- 
tors sampled and reduces the statistical noise by a large 
factor for the case of pseudoscalar mesons. We also take 
a random starting point for our time sources for the very 
coarse, coarse and fine ensembles. 

The pseudoscalar meson correlation function C a b(t) for 
meson of valence content ab is calculated by multiplying 
together the quark propagator for quark a and the com- 
plex conjugate of the quark propagator for quark b from 
the same source on a given configuration, matching col- 
ors at the source and sink and matching the sink spatial 
site index, which is summed over to set the meson to zero 
momentum. The meson correlation function is then av- 
eraged over time sources on a single configuration. This 
means that any correlations between the time sources on 
a given configuration are accounted for. We also have 
to worry about autocorrelations between results on suc- 
cessive configurations in an ensemble. Tests by binning 
correlators have shown that the results on different con- 
figurations are independent of each other except on the 
finest lattices. We therefore bin the correlators on su- 
perfine and ultrafine lattices by a factor of two. 

The correlation function averaged over the indepen- 
dent samples from an ensemble is then fit as a function 
of the time separation between source and sink, t, to the 
form: 



C(t) 



We" 



-Mit 



-Mt(T-th 



(8) 



for the case a = b. i = is the ground state and larger 
i values denote radial or other excitations with the same 
J quantum numbers. T is the time extent of the lat- 
tice. For the unequal mass case there are additional 'os- 
cillating' terms coming from opposite parity states, de- 



5 



noted i p : 

C(t) = ^e' Aht + (-lYa ip e- M ^ + (t->T-t) (9) 

To fit we use a number of exponentials i, and where 
appropriate i p , in the range 2-6, loosely constraining the 
higher order exponentials by the use of Bayesian pri- 
ors [33]. As the number of exponentials increases, we 
see the \ 2 value fall below 1 and the results for the fit- 
ted values and their errors for the parameters for the 
ground state i — stabilise. This allows us to determine 
the ground state parameters ao and Mq as accurately 
as possible whilst allowing the full systematic error from 
the presence of higher excitations in the correlation func- 
tion. We take the fit parameters to be the logarithm of 
the ground state masses Mq and Mo and the logarithms 
of the differences in mass between successive radial ex- 
citations (which are then forced to be positive). The 
Bayesian prior value for Mo is obtained from a simple 
'effective mass' in the correlator and the prior width on 
the value is taken as a factor of 1.5. The prior value for 
the mass splitting between higher excitations is taken as 
roughly 600 MeV with a width of a factor of 2. Where 
oscillating states appear in the fit, the prior value for 
Mq p is taken as roughly 600 MeV above Mo with a prior 
width of a factor of 2 and the splitting between higher 
oscillating excitations is taken to be the same as for the 
non-oscillating states. The amplitudes ai and at are 
given prior widths of 1.0. 

Our fit includes the effect of correlation between dif- 
ferent values of t. We apply a cut on the range of eigen- 
values from the correlation matrix that are used in the 
fit of 10~ 3 or 10~ 4 . We also cut out very small t values 
from our fit, typically below 3 or 4, to reduce the effect 
of higher excitations. 

The results for masses and amplitudes from fits in 
equations [8] and [9] are in units of the lattice spacing. The 
value of the lattice spacing must be determined for each 
ensemble to enable conversion to physical units. For this 
we use the parameter r-y, defined from the heavy quark 
potential [4|. This parameter can be accurately deter- 
mined (to better than 0.5%) in units of the lattice spacing 
and so is good for making ensemble to ensemble compar- 
isons of a. Results from the MILC collaboration are given 
in Table [TJ Unfortunately n does not have a directly ac- 
cessible physical value. Instead we must determine that 
from other quantities. In [31] we used four other physical 
quantities with well-known experimental values to fix the 
value of 7"i to 0.3133(23) fm. This then yields a value for 
a on a given ensemble with two errors - an error from the 
value of T\ja on that ensemble and an error, correlated 
between ensembles, from the physical value of T\. 

The amplitude, ao, from the fits in equations [8] and [9] 
is directly related to the matrix element for the local 
pseudoscalar operator to create or destroy the ground- 
state pseudoscalar meson from the vacuum. From the 
definition of the correlator and using a relativistic nor- 



malisation for the fields: 

a = ((Q\P S \P)) 2 /2M (10) 

where the pseudoscalar current P s — 0756 for quark con- 
tent ab. Because of the chiral symmetry of the staggered 
quark action we have a partially conserved axial current 
(PCAC) relation that relates the local pseudoscalar oper- 
ator above to a temporal axial current that is absolutely 
normalised on the lattice. This allows us to determine 
the decay constant for these pseudoscalar mesons with- 
out having to worry about an uncertainty from the renor- 
malisation between the lattice and the continuum. The 
decay constant for meson P with quark content ab is de- 
fined from: 

(Q\a^ l5 b\P(p)) = f PPfi . (11) 

For a meson at zero momentum, and using the PCAC 
relation d^A^ = (rn a + nib)P s to relate the axial vector 
and pseudoscalar currents, this becomes: 

(m, + m b )(0\a l5 b\P(p)) = f P M 2 P , (12) 

where m a and nib are the appropriate quark masses. 
Combining this with equation [10] then allows us to deter- 
mine fp in lattice QCD from our fits to the correlators 
for pseudoscalar meson P using 



fp = (m a + m b )y ^3. (13) 

Here m a and nib are the quark masses used in the lattice 
QCD calculation. 

fp in turn is related, for charged pseudoscalars such 
as the 7r, K, D and D s mesons, to the experimentally 
measurable leptonic branching fraction via a W boson: 

B(P^lM"f)) = ' fpmfmpll i- , 

air \ m P ) 

(14) 

up to calculable electromagnetic corrections. V ab is the 
appropriate CKM element for quark content ab. Tp is 
the pseudoscalar meson lifetime. 

III. RESULTS 

Accurate results for the D s meson require accurate 
tuning of both the c and the s quark masses. We use 
the pseudoscalar mesons made purely of c quarks or of 
s quarks to do this and so first discuss results for these 
mesons. 

Table [TT] lists the valence HISQ quark masses close to 
that of the charm quark that we used for each of the 
gluon configuration ensembles along with the correspond- 
ing value of the Naik parameter (1 + e). We also list the 
values of the ground-state pseudoscalar cc meson mass 
and decay constant obtained from our fits to the cc me- 
son correlators to equation [8] The decay constant, f Ve , 
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will be discussed in subsection III C - it is a useful quan- 



tity to calculate despite the fact that the r/ c is a neutral 
particle and does not undergo a purely leptonic decay of 
the kind given in equation |14| To tune the charm quark 
mass in the HISQ action we must interpolate to the point 
at which the mass of the r\ c has the correct physical value 
on each ensemble. This physical value is not exactly the 
experimental value (2.980 GeV [35]) because our lattice 
QCD calculation corresponds to a world in which there 
are no electromagnetic interactions and we do not allow 
our i] c meson to annihilate to gluons. Both of these ef- 
fects tend to reduce the rj c mass by small amounts and 
so the appropriate physical value for us to compare our 
lattice QCD calculation to is 2.985(3) GeV, allowing a 
50% error for each correction to the experimental value. 
The corrections are obtained from a potential model for 
the electromagnetic effect and from perturbation theory 
for the effect of gluon annihilation [34] . 

Figure [T] shows the meson mass in physical units plot- 
ted against the quark mass, also in physical units, for 
each ensemble. This plot demonstrates how the quark 
mass tuning is done, as well as illustrating very clearly 
how accurately we can do this from lattice QCD. Several 
features of the figure stand out. On a given ensemble 
the value of the meson mass is linear in the quark mass, 
as we would expect. The lines showing this behaviour 
(not plotted on the figure) are essentially parallel with 
a slope close to the naive expectation of 2 for ensembles 
with different lattice spacing values. In fact the slope 
does increase from 1.7 on the very coarse lattices to 2.3 
on the superfine lattices. The reason for this is that the 
x-axis is a well-defined 'running' quark mass, being the 
quark mass in the HISQ Lagrangian with a particular ul- 
traviolet scale set by the lattice spacing. This is why we 
denote the mass on the x-axis as m c (a). The horizontal 
line indicates the correct value of the r\ c mass and there- 
fore, where it cuts each set of results, the tuned value 
of m c at that lattice spacing. These values 'run' to the 
left on finer lattices as the ultraviolet cut-off increases, 
as expected from perturbation theory. We expect the 
variation of r\ c mass with quark mass to be some number 
(say, 2) times the quark mass at a fixed scale. Therefore 
on finer lattices, where the scale is higher, we expect the 
slope to be larger, as demonstrated in Figure [TJ 

Another feature is that the results for different ensem- 
bles with very similar values of the lattice spacing are 
very close together i.e. there is very little dependence of 
the tuned c mass on the sea quark masses. The results 
for different physical volumes (sets 4 and 5) lie on top of 
each other showing that there is no dependence on the 
volume. We would not expect any significant volume de- 
pendence on these large spatial volumes for the r\ c since 
it is a relatively small particle. 

From the horizontal line on Figure [T] and the lattice 
points on the line it is clear that we have tuned the charm 
quark mass very well on all except the superfine lattices 
(where it is off by 0.1%). In each case this corresponds to 
the lightest charm quark mass in our Table [IT] Figure [I] 



> 
o 




0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 
m c (a) (GeV) 



FIG. 1: Results for the mass of the pseudoscalar meson made 
of quarks with masses close to that of the charm quark mass 
for the full set of ensembles from Table [I] The a>axis is the 
lattice bare mass of the quark, which runs with lattice spacing 
from right to left. Very coarse ensembles are triangles; coarse, 
crosses; fine, squares; superfine, pluses; ultrafine, bursts. Re- 
sults for heavier sea u/d quark masses at each lattice spacing 
are in red, lighter ones are in blue. On the coarse lattices the 
very heavy sea masses of set 6 are in pink, the lighter strange 
sea mass of set 7 in grey and the large volume results on set 5 
are in green, on top of the result from set 4. Statistical errors 
are too small to be visible on this plot. The results show that 
tuning the quark mass to that of charm depends very little 
on the sea quark masses or on the volume. The dotted line 
gives the physical value, with its error, appropriate to lattice 
QCD, see text. 



does not include errors in converting the lattice quark 
mass or rj c mass to GeV coming from the values of r±/a 
or the physical value of r\. The effect of these errors is 
reduced over naive expectations because Am ric is close 
to 2Ato c , and so the leading order change from any Aa 
cancels out. This issue was addressed in |36j . Here we 
are not aiming to determine m c , but simply to make sure 
we understand the errors in other quantities induced by 
the tuning error in m c a, so we leave a more detailed dis- 
cussion of this source of systematic error to the sections 
on the individual quantities. 

Table [TT| lists the valence HISQ quark masses close to 
that of the strange quark that we used for making strange 
quark propagators on each of the gluon field ensembles. 
We also list the corresponding values of the mass of the 
ground-state ss meson, the rj s . The r] s is not a particle 
available to study in the real world where ss mixes with 
uu and dd. However, by omitting these possibilities in the 
lattice QCD calculation we can obtain a pure ss 'pion- 
like' meson. This turns out to be useful for tuning the 
s quark mass because the r] s mass can be determined 
relatively precisely, and is less sensitive to the sea quark 
masses than, for example, m v . However, the physical 
value for the rj s mass has to be determined by relating 
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Set 


am a 


1 + e 


om, c 




am a 


am Vs 


1 


0.81 


0.665 


2.19381(16) 


0.3491(5) 


0.061 


0.50490(36) 




0.825 


0.656 


2.22013(15) 


0.3539(5) 


0.066 


0.52524(36) 




0.85 


0.641 


2.26352(15) 


0.3622(5) 


0.080 


0.57828(34) 


2 


0.825 


0.656 


2.21954(13) 


0.3537(4) 


0.066 


0.52458(35) 


3 


0.622 


0.779 


1.79132(8) 


0.25706(18) 


0.0489 


0.41133(17) 




0.65 


0.762 


1.84578(8) 


0.26368(18) 


0.0537 


0.43118(18) 


4 


0.63 


0.774 


1.80849(11) 


0.25998(20) 


0.0492 


0.41436(23) 




0.66 


0.756 


1.86666(10) 


0.26721(20) 


0.0546 


0.43654(24) 




0.72 


0.720 


1.98109(10) 


0.28228(22) 


0.06 


0.45787(23) 




0.753 


0.700 


2.04293(10) 


0.29114(24) 


0.063 


0.46937(24) 


5 


0.63 


0.774 


1.80856(7) 


0.26006(15) 


0.0492 


0.41457(14) 


6 


0.625 


0.777 


1.79347(13) 


0.2556(3) 


0.0491 


0.41196(24) 












0.0525 


0.42588(30) 












0.0556 


0.43834(30) 


7 


0.619 


0.781 


1.78595(15) 


0.2564(3) 


0.0487 


0.41030(31) 


8 


0.413 


0.893 


1.28057(7) 


0.17217(11) 


0.0337 


0.29413(12) 




0.43 


0.885 


1.31691(7) 


0.17508(11) 


0.0358 


0.30332(12) 




0.44 


0.880 


1.33816(7) 


0.17678(11) 


0.0366 


0.30675(12) 




0.45 


0.875 


1.35934(7) 


0.17850(11) 


0.0382 


0.31362(14) 


9 


0.412 


0.894 


1.27522(7) 


0.17086(10) 


0.0336 


0.29309(13) 




0.427 


0.885 


1.30731(10) 


0.17344(15) 


0.03635 


0.30513(20) 


10 


0.273 


0.951 


89935('12'l 


1 1864l'24 > l 


0228 


20621 (IQl 




0.28 


0.949 


0.91543(8) 


0.11986(21) 


0.024 


0.21196(13) 


11 


0.193 


0.975 


0.66628(13) 


0.0882(3) 


0.0161 


0.15278(28) 




0.195 


0.975 


0.67117(6) 


0.08846(11) 


0.0165 


0.15484(14) 












0.018 


0.16209(17) 



TABLE II: Results for the masses in lattice units of the goldstone pseudoscalars made from valence HISQ charm or strange 
quarks on the different MILC ensembles, enumerated in Table [I] Columns 2 and 3 give the corresponding bare charm quark 
mass, and Naik coefficient respectively. Column 6 gives the bare strange quark mass (e = in that case). A lot of the meson 
masses in this table appear also in [55] but we have added results on the coarse 02/05 and 01/03 ensembles (sets 6 and 7) and 
the large volume coarse 01/05 ensemble (set 5) as well as improving the tuning of masses on other ensembles and improving 
some fits on sets 4, 10 and 11. Results for the decay constant of the r) c meson are also included, for analysis in subsection |III C| 



it to 7r and K meson masses known from experiment. 
In earlier lattice QCD calculations we determined m Va 
= 0.6858(40) GeV 134] and this is the value we will use 
here. We also studied the rj s decay constant which is 
again a quantity that cannot be measured experimentally 
in the real world but one which turns out to be useful for 
determining the lattice spacing. We will not discuss f Vg 
further here. 

Figure [2] shows the square of the r] s mass against the 
quark mass, both in physical units, for each ensemble. 
We expect oc m s from leading order chiral perturba- 
tion theory and the results indeed show this dependence. 
Once again the lines demonstrating this (not plotted on 
the figure) are fairly parallel but with a slope increasing 
on the finer lattices as the quark mass for a given me- 
son mass runs to smaller values. The horizontal plots 
gives the physical value of the rj s mass given above and 
the strange quark mass can be read off for each ensem- 
ble from where this crosses the line of data. Again we 
have well-tuned strange quark masses at each value of 
the lattice spacing at the lightest end of the range. The 
strange quark mass values on the very coarse and coarse 
lattices are rather close together but on the finer lat- 
tices the strange mass changes as rapidly with lattice 



spacing as the charm mass does in Figure [T] In the con- 
tinuum limit the ratio of these two masses becomes a 
scale-invariant constant [36] . 

Again it is evident from Figure [2] that there is very lit- 
tle dependence of the tuned s quark mass on either the 
sea quark masses or the volume. Because the value of 
the tuned s quark mass is proportional to the square of 
the r] s mass the relative uncertainty in m s arising from 
lattice spacing errors is equal to that of the lattice spac- 
ing. There is no cancellation as there was in the case of 
the charm quark. In addition the 0.6% uncertainty in the 
physical value of the r) s mass is significant, because it be- 
comes an uncertainty of 1.2% in m s . The effect of these 
uncertainties on the mass and decay constant of the D s 
meson will be discussed below. 

The staggered quarks in the sea are asqtad improved 
staggered quarks rather than HISQ quarks, i.e. they use 
a different discretisation of the quark piece of the QCD 
Lagrangian. The s quark mass in the two formalisms 
will then not be the same, and we need to understand the 
ratio of the two so that we can extrapolate to the physical 
(real world) point for both the valence and sea quark 
masses. We can determine the physical points for the sea 
quark masses from our tuning of the valence masses and 
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FIG. 2: Results for the square of the mass of the pseudoscalar 
meson made of quarks with masses close to that of the strange 
quark mass for the full set of ensembles from Table [I] Errors 
are statistical errors from the fits to the meson correlators. 
The a;-axis is the lattice bare mass of the quark, which runs 
with lattice spacing from right to left. Very coarse ensembles 
are triangles; coarse, crosses; fine, squares; superfine, pluses; 
ultrafine, bursts. Results for heavier sea u/d quark masses at 
each lattice spacing are in red, lighter ones are in blue. On 
the coarse lattices the very heavy sea masses of set 6 are in 
pink, the lighter strange sea mass of set 7 in grey and the 
large volume results on set 5 are in green, on top of the result 
from set 4. The results show that tuning the quark mass to 
that of strange depends very little on the sea quark masses or 
on the volume. The dotted line gives the physical value, with 
its error, appropriate to lattice QCD, see text. 



this ratio. There is very little sea quark mass dependence 
in the quantities that we study here, so that we do not 
need to know this ratio accurately. It is discussed further 
in Appendix |A"} 

Once we have determined the c and s masses to be used 
to give the required physical results for the r\ c and the r\ s 
mesons, the D s meson correlator is entirely prescribed. 
There are no further adjustable parameters, given the na- 
ture of QCD. The fit to the D s meson correlators gives us 
both the D B meson mass (from M in equation [9]) and its 
decay constant (from clq) as testable outputs from lattice 
QCD. Since the D s meson mass is well-known experi- 
mentally it provides an excellent independent test of the 
error analysis on the decay constant. It is therefore very 
important to analyse both of these quantities together. 



The D s meson correlators are made from the same c 
and s quark propagators that are used for the r\ c and 
rj s above. We must use equation [9] to fit the D s correla- 
tors, however, because they do have additional oscillating 
terms in them. Table IIIII lists results for the masses, Mq 



TABLE III: Results for the mass and decay constant of the 
D s meson in units of the lattice spacing for a range of charm 
and strange quark masses on each MILC ensemble. 



Set 


am c 


am s 




afo s 


1 


0.81 


0.061 


1.4665(8) 


0.1970(10) 




0.825 


0.066 


1.4869(7) 


0.1994(10) 




0.825 


0.080 


1.5019(6) 


0.2042(8) 




0.85 


0.066 


1.5117(8) 


0.2004(10) 




0.85 


0.080 


1.5266(6) 


0.2053(9) 


2 


0.825 


0.066 


1.4869(11) 


0.1997(20) 


3 


0.622 


0.0489 


1.1890(7) 


0.1538(9) 




0.65 


0.0537 


1.2247(5) 


0.1561(9) 


4 


0.63 


0.0492 


1.2007(5) 


0.1559(7) 




0.66 


0.0546 


1.2391(5) 


0.1586(6) 




0.66 


0.06 


1.2452(5) 


0.1604(6) 




0.66 


0.063 


1.2486(4) 


0.1614(6) 




0.72 


0.0546 


1.3027(6) 


0.1602(7) 




0.72 


0.06 


1.3086(5) 


0.1620(7) 




0.72 


0.063 


1.3120(5) 


0.1631(6) 




0.753 


0.0546 


1.3369(6) 


0.1610(7) 




0.753 


0.06 


1.3429(5) 


0.1629(7) 




0.753 


0.063 


1.3462(5) 


0.1639(7) 


5 


0.63 


0.0492 


1.2013(5) 


0.1561(8) 


6 


0.625 


0.0491 


1.1916(7) 


0.1553(10) 


7 


0.619 


0.0487 


1.1867(10) 


0.1548(17) 


8 


0.413 


0.0337 


0.84721(23) 


0.10836(24) 




0.43 


0.0358 


0.86982(23) 


0.10943(24) 




0.43 


0.0366 


0.87079(22) 


0.10970(24) 




0.43 


0.0382 


0.87274(21) 


0.11028(24) 




0.44 


0.0358 


0.88152(23) 


0.10959(27) 




0.44 


0.0366 


0.88249(23) 


0.10986(27) 




0.44 


0382 


88443(221 


0.11044(24) 




0.45 


0.0358 


0.89317(24) 


0.10974(27) 




0.45 


0.0366 


0.89414(23) 


0.11001(27) 




0.45 


0.0382 


0.89607(23) 


0.11059(27) 


9 


0.412 


0.0336 


0.84352(26) 


0.10779(31) 




0.427 


0.03635 


0.86443(40) 


0.1086(5) 


10 


0.273 


0.0228 


0.59350(24) 


0.07500(27) 


11 


0.193 


0.0161 


0.43942(33) 


0.05533(39) 




0.195 


0.0165 


0.44270(28) 


0.05550(34) 



and the decay constant derived from ao for each com- 
bination of c and s masses that we have used on each 
ensemble. The statistical errors coming from the fit are 
significantly larger for the D s than for the rj c . This is 
because the noise in heavy-light correlators has a lower 
mass associated with it than the signal. The mass in the 
squared correlator which gives the noise is given by one 
half of the sum of the r\ c and r\ s masses, which is smaller 
than the signal D s mass. This means that the signal to 
noise ratio degrades at large times for the D s correlator 
and the statistical error increases. This is illustrated in 
Figure [3] in which we explicitly plot and compare the 'ef- 
fective mass' extracted from the D s correlator and from 
its statistical error. This issue becomes a problem for B 
meson correlators [37] ■ It is not a big problem for the 
D s , however, and the statistical errors that we obtain in 
Table [nT] are very small. 

To determine the physical mass of the D s meson as ac- 
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FIG. 3: Results for the effective mass of the D a correlator and 
the effective mass of the noise in the D a correlator plotted as 
a function of lattice time for one correlator on the fine lattices 
(set 8). The effective mass is obtained from the log of the ratio 
of the correlator (or its error) at successive times. At large 
times it becomes the mass of the lowest state in the correlator 
or its error. The lines compare the results to the expected 
mass i.e. the D a mass for the signal and (m Vs + m Vc )/2 for 
the noise. 



curately as possible we want to minimise errors coming 
from the conversion from lattice units to physical units 
i.e. from the lattice spacing. The error on the physical 
value of 7"i is 0.7%. Applied directly to the D s mass this 
would amount to a sizeable 14 MeV error. This can be 
avoided however, by calculating instead the mass differ- 
ence itid b — m r)e /2. Because this is much smaller (480 
MeV) it will have a much reduced absolute error from 
the lattice spacing [3J. In addition, it is much less sen- 
sitive to any errors from mistuning of the c quark mass 
because the leading contribution of m c effectively cancels 
in this difference. Indeed this difference can be thought 
of as the difference in binding energy between a charmo- 
nium meson and a charm-light meson, and is therefore 
an important physical quantity. The fact that it can be 
calculated accurately in lattice QCD and compared to 
experiment is a stringent test of QCD itself. 

The first stage in the analysis of the D s meson mass 



is to determine the difference mo s 



„ /2 for tuned 



c and s quark masses on each ensemble. As discussed 
above, we have results very close to the tuned point on 
almost every one of the 11 ensembles. However, it is 
important to make sure that all of our results are tuned 
to the same point before extrapolation and so we first 
test the dependence of mj, — m, )c /2 as a function of m s 
and m c . Figures E] and [H] show results as a function of 
and m Vc for sets 1, 4 and 7 where we have multiple data 
points with different combinations of m c and m s and so 
can unravel the separate dependences. The dependence is 
plotted against meson mass rather than directly against 
the quark mass since the tuning condition is set by the 
rj c or rj s meson mass, so this is a more direct (and more 
physical) way to study any mistuning effects. Note that 
the mass values of the rj c and rj s are above their physical 
values for the cases given in figures [4] and [5] Since we are 
only studying small mistuning effects for the values of the 



0.54 



0.53 



| 0.52 



r , 0.51 - 



0.5 



0.49 









i i — 


■••;>-''' 
































11 
X 


•X 

,a'S-" 
a' , 
.11 


..in"" 
Hi 
hi 


■ ■ 





0.48 0.5 0.52 0.54 0.56 0.58 0.6 0.62 0.64 
m^ 2 (GeV 2 ) 



FIG. 4: Results for the mass of the D 3 meson (specifically 
the difference between that mass and one half of the rj c mass) 
as a function of the square of the rj a meson mass, acting as a 
proxy for the strange quark mass. Results are for a range of 
different quark masses around the masses of the c and s quark 
masses on very coarse set 1 (triangles), coarse set 4 (crosses) 
and fine set 7 (squares). The lines are fits to the results for 
each ensemble allowing linear terms in m^ s and m IJc . Here 
the lines join points for a fixed c quark mass. See Figure[5]for 
the equivalent as a function of m Vc . 
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FIG. 5: Results for the mass of the D s meson (specifically 
the difference between that mass and one half of the rj c mass) 
as a function of the r\ c meson mass, acting as a proxy for 
the charm quark mass. Results are for a range of different 
quark masses around the masses of the c and s quark masses 
on very coarse set 1 (triangles), coarse set 4 (crosses) and 
fine set 7 (squares). The lines are fits to the results for each 
ensemble allowing linear terms in m^ s and rn rla . Here the 
lines join points for a fixed s quark mass. See Figure [4] for the 
equivalent as a function of m^ s . 
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masses that we have closer to the physical points, this will 
give a sufficiently accurate picture of these effects. 

In figure [4] we see that the dependence of tud, — m ri j2 
on is linear as we expect, since this corresponds to 
a linear dependence on m s . The slope is clearly physi- 
cal i.e. independent of the lattice spacing (whereas the 
slope against m s would not be, because of the running 
of m s itself, discussed earlier). The value of the slope is 
0.20(1) and this can be compared to an 'experimental' 
slope, albeit over a much larger mass range, of 0.22 ob- 
tained by comparing results for the masses of the D and 
the D s [35] . Figure [5] also shows linear dependence on 
m c , expressed physically as linear dependence on m Vc . 
The slope does differ on the very coarse lattices from the 
others so showing some lattice spacing dependence in this 
case. The slope is also very small ~ 0.05 because, as dis- 
cussed above, the leading dependence on m c cancels be- 
tween rriD a and m^/2. The slope is again similar to the 
'experimental' value of 0.03 obtained over a much larger 
mass range from comparing B s and D s mesons |35j . 

Results from Figures [4] and [5] can be used to adjust the 
values of to_d s — m Vc /2 on each ensemble to the tuned 
point, m Vc = 2.985 "GeV and m Vs = 0.6858 GeV. An 
error of 50% of any shift is added in quadrature to the 
statistical error. The shifts from mistuning are less than 
the statistical error on all ensembles except sets 2 (very 
coarse) and 10 (superfine). On set 10 the shift is by 
1.5 times the statistical error and on set 2 by 4 times 
the statistical error. Table |IV| gives the tuned value of 
m]j 3 — m ric /2 in GeV on each ensemble along with two 
errors. The first is the statistical/tuning error and the 
second is that from the error in r±/a on that ensemble. 
This error is a factor of 3 smaller than its naive value 
because of a cancellation of lattice spacing errors inside 
the mass difference. Any change in r\/a means a change 
to m Vc and m Vs as well as a change in m£) s — m Vc /2. The 
results then need to be retuned to the physical c and s 
masses and this largely cancels the change resulting from 
the change in r±/a. The error from r\/a uncertainty is 
much smaller than the statistical error then in every case. 
The statistical errors, which dominate, are at the level of 
1 MeV. 

We can then extrapolate the tuned values on each en- 
semble in the lattice spacing and the sea quark masses to 
the physical point where the lattice spacing is zero and 
the sea quark masses take their real world values. It is 
clear from Table |TV] looking at the coarse and fine ensem- 
bles that m,D s — va r)a j2 has no significant dependence on 
the sea quark masses at the level of our 1 MeV statistical 
errors. The picture is obscured on the very coarse lattices 
by the larger error on set 2 from mistuning. In fact if we 
compare sets 1 and 2 at the rj c and rj s masses correspond- 
ing to those available on set 2 (i.e. at somewhat heavier 
masses than the correctly tuned point) then we find again 
that sets 1 and 2 agree on the value of mo 3 — m ric but 
now within an error of 1.5 MeV rather than the 3 MeV 
in Table [TV] 

We expect mo 3 —m nc /2 to be very insensitive to the sea 



TABLE IV: Values for the mass and decay constant of the 
D 3 meson and for the decay constant of the rj c after tuning 
to the physical c and s masses (i.e. the physical r\ c and rj 3 
meson masses) on each ensemble. Results are in GeV with 
two errors, the first from statistics and tuning and the second 
from the uncertainty in r\/a on that ensemble. 



Set 


m Ds _ (GeV) 


fn s (GeV) 


L (GeV) 

J r ic V / 


1 

2 


0.5021(12)(4) 
0.5020(32) (4) 


0.2674(14) (3) 
0.2671(28)(3) 


0.4753(9)(2) 
0.4756(6)(2) 


3 
4 
5 
6 
7 


0.4889(12)(3) 
0.4897(9) (3) 
0.4906(9) (3) 
0.4909(12) (3) 
0.4911(17)(3) 


0.2564(15)(2) 
0.2573(12)(2) 
0.2576(13)(2) 
0.2586(17)(2) 
0.2592(28)(2) 


0.4284(3) (1) 
0.4291(4)(1) 
0.4292(4)(1) 
0.4255(5)(1) 
0.4286(6)(1) 


8 
9 


0.4823(6) (2) 
0.4817(6)(2) 


0.2525(6)(2) 
0.2520(7)(2) 


0.4012(3)(2) 
0.3998(3)(2) 


10 


0.4784(10) (2) 


0.2499(9)(3) 


0.3945(10)(3) 
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0.4766(13) (4) 


0.2481(17)(5) 


0.3953(13)(6) 



quark masses based on chiral perturbation theory. This 
couples a nonrelativistic lagrangian for D and D s meson 
fields to the pion octet and gives an expansion in powers 
of 7r, K and ?78 masses for the mass and decay constant 
of the appropriate D meson. The D s has valence c and s 
masses which have been tuned to the appropriate values 
so the only dependence we are interested in here is the 
dependence on sea s and u/d quark masses which enter 
through the masses of mesons made either purely of sea 
quarks or of mixed sea and valence quarks. The lead- 
ing tree-level dependence on sea quark masses is a term 
C(2mi tSea +m SjSea ). Loops couple the D s meson to a vir- 
tual DK or D s rj s pair. This generates logarithmic terms 
but with, in this case, a very benign dependence on sea 
quark masses since none of the associated meson masses 
vanish in the chiral limit. These terms can then simply 
be viewed as additional polynomial terms in mi tSea and 
m SjSea . A more detailed chiral analysis is not useful here 
because the sea quark mass dependence of our results is 
clearly so small as to have no useful information in it. We 
simply need to make sure that we allow a sufficient error 
on the extrapolated value at the physical point to allow 
for any sea quark mass dependence that might be there. 
For this purpose a simple polynomial expansion in mi tSea 
and m S)Seo suffices. We take as expansion coefficients Sxi 

and Sx s where SXq — (^Tlg.sea TTlq,sea,phys^)/'ms,sea,phys- 

ms,sea,phys is the sea (asqtad) strange quark mass at the 
physical point. We take this value from results quoted by 
the MILC collaboration [27 for very coarse to superfine 
and use the analysis of the ratio of HISQ to asqtad mases 
from Appendix [A] to give the value of rn s ^ sea:P hys on the 
ultrafine lattices. We take m Lsea , phys = m s ^ ea ^ hys /27.2 
using the ratio for mi/m s determined by the MILC col- 
laboration [57] ■ Table |Tv| shows that m,D s — m nc /2 does 
have significant dependence on the lattice spacing, chang- 
ing by 20 times the statistical error between very coarse 
and ultrafine lattices. This is also not surprising because 
the charm quark is relatively heavy and consequently the 
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scale for discretisation errors here will be much higher 
than that for quantities involving only light quarks. This 
is why it is important to have a formalism, such as HISQ, 
with very well controlled discretisation errors and to have 
results at many values of the lattice spacing. Discretisa- 
tion errors with the HISQ action can appear only as pow- 
ers of a 2 - no odd powers of a are allowed. The a 2 errors 
appearing at tree level have been removed and so the co- 
efficient of a 2 terms is O(ot s ). The inclusion of the Naik 
term with coefficient calculated at tree-level means that 
all (m c a) 2n discretisation errors are removed at leading 
order in v 2 /c 2 where v 2 is the velocity of the charm quark 
in the D s or rj c . Thus discretisation errors from the HISQ 
action are expected to be at the level of (v 2 /c 2 )(m c a) 2n , 
except for the a 2 term which is further suppressed by a s . 
There are additional a s a 2 and tree-level a 4 and higher 
errors coming from the gluon action, however. These we 
would typically expect to have a scale of a few hundred 
MeV (i.e. Aqcd) associated with them rather than m c , 
so their effects will be included if we allow for a scale of 
m c . 

We therefore take the following fit form to extrapolate 
A = m]j 3 — m rk /2 to the physical point: 



4 

A(a,Sxi,Sx s ) = Ap h y S [1 + ^ cj (m c a) 



) 2 :i 



(15) 



+ 2bidxi(l + c b (m c a) 2 ) 
+ 2b s 5x s (l + c s (m c a) 2 ) 
+ 4bu(5xi) 2 +2bi s 5xi6x s 



b sa (5x s y 



We use a constrained fit [33J to this form which allows 
us to estimate the errors arising from different pieces of 
the fit. The prior value and width for A p h ys we take 
as 0.5, with the very broad width of 0.2. Note that we 
give the discretisation errors a scale of m c . The prior 
value and width that we take on the c„ parameters is 
0.0(2), estimating v 2 /c 2 for the c quark inside the D s to 
be 0.2. Ci, which multiplies the a 2 errors, is a factor of ot a 
smaller from the arguments above so we take the prior 
for ci to be 0.00(6). The b parameters multiplying the 
linear sea quark mass dependence are taken to have prior 
values and widths of 0.00(7). The size of the prior width 
here is set by the fact that the dependence of A on the 
valence light quark mass inside the D s is known from a 
comparison of D and D s . This would give a slope with 
valence mass, in units of the strange mass, of 0.2. Sea 
quark mass effects are a factor of at least 3 smaller than 
valence mass effects in gold-plated quantities, so we take 
a prior width of 0.07. By the same reasoning we allow 
the b parameters multiplying the quadratic dependence 
to be as large as (0.2) 2 /3, i.e. we take the prior on these 
parameters to be 0.000(13). 

The extrapolated result at the physical point, A p h ys 
from the fit above is 0.4753(22) MeV with a x 2 /dof of 
0.2 for 11 degrees of freedom. We fit all of the data 
including the two volumes for the coarse lattices, sets 4 
and 5. Missing out set 5 makes no appreciable difference 



to the result. Modifications to the fit form above also 
do not change this number significantly. Here we itemize 
the effect of some of them: 

• changing the prior on all Cj (including ci) to 0.0(5) 
changes A p h ys 0.4<r and increases the error by 20 

• adding two extra powers of a 2 into the sum on j in 
equation 15 (i.e. using 6 terms instead of 4) does 
not change A p h ys or the error at all. The same is 
true for subtracting two powers of a 2 (i.e. using 2 
terms instead of 4). 

• adding extra discretisation errors into the sea-quark 
mass dependence (i.e. a term proportional to 
(m c a) 4 in each of the terms linear in 8xi and 5x s 
and a term proportional to (m c a) 2 in each of the 
quadratic terms) makes no difference at all. 

• missing out the sea quark mass dependence alto- 
gether makes no difference to A p h ys but increases 
the x 2 value to 0.33. 

• Changing all the 8x values by 10% in either di- 
rection makes no appreciable difference, nor does 
changing them within their error bars on, for ex- 
ample, the ultrafine or fine lattices. 

• missing out the very coarse lattice results makes 
no difference; missing out the very coarse and the 
coarse shifts A p h ys by 0.4<r (1 MeV), and increases 
the error to 3 MeV as x 2 drops to 0.1. 

• missing out the ultrafine result shifts A p ^j, s by 0.2ct 
(0.5 MeV) and increases the error to 3 MeV. 

Figure [6] shows the results plotted against the square 
of the lattice spacing along with the fitted curve above, 
taken at the physical sea quark mass values (i.e. Sxi = 
Sx s = 0). The value plotted on the y-axis is rriD s itself, 
generated by adding m, la /2 = 1.4925 GeV to A. The 
result at a = is then the value of the D s mass in a 
world without electromagnetism. To compare to exper- 
iment we need to estimate and add in the effect of the 
electromagnetic repulsion between the positively charged 
quark and antiquark inside the D s . To do this we com- 
pare experimental masses for the D + , D°, D s , B + , B° 
and B s to a phenomenological formula allowing for elec- 
tromagnetic effects proportional to the product of quark 
and antiquark electric charges inside the meson as well 
as the square of the electromagnetic charge on the light 
quark. This latter term is a self-energy effect, not needed 
for the heavy quarks because it will cancel in all the dif- 
ferences taken (and therefore is absorbed into the heavy 
quark mass). In comparing charged and neutral mesons 
containing u and d quarks we must allow for the mass dif- 
ference between u and d quarks. Then we can write [38j : 

M(Q, q) = M slm {Q 1 q)+Ae q e Q +Be 2 q +C{m q -m l ) (16) 

where M S i m is the mass of the meson in the absence of 
electromagnetism and with m u = m&. If we take exper- 
imental results for the meson masses above along with 
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FIG. 6: Results for the mass of the D s meson tuned to the 
correct valence c and s mass on each ensemble from Table [iVl 
as a function of the square of the lattice spacing. The line 
shows the result of the fit described by equation |15| taken at 
the physical values for the sea quark masses. The shaded band 
gives our final result adjusted for electromagnetic effects and 
with the full error as described in the text. The black burst 
gives the experimental result. 

Tn a /mi — 27.2 and m u /m& — 0.42 we obtain A « 4 
MeV, B w 3 MeV and Cm s m 100 MeV. The latter 
quantity differs by 10% between D and B mesons, indi- 
cating 1/mg effects at this level that we ignore here. The 
resulting electromagnetic shift for the D s is then 1.3(7) 
MeV, where we take an error of 50% on the shift, safely 
encompassing 1 /mg effects and other limitations of this 
model. Adding 1.3 MeV to our fit result gives the shaded 
band in Figure [6j where we now include our full error of 
3.2 MeV. The full error budget is discussed below. 

Figure [7] shows the sea quark mass dependence of our 
results plotted against 5x[. The fitted curves are those 
from equation [l5j For each group of ensembles we use 
the lattice spacing value from the ensemble with light- 
est sea quark mass to plot the fit curve. No significant 
dependence on Sxi or 5x s is evident. 

Table |V| shows the complete error budget for mu s from 
our calculation. The error of 2.2 MeV from our fit to A 
above includes the effect of statistical errors (including 
valence mass mistuning errors), T\/a errors and errors 
arising from the extrapolation in sea quark masses and 
lattice spacing. We can separate these errors as described 
in 31J by working out how the final error changes when 
any of the inputs to the fit changes and dividing a 1 into 
a sum of terms coming from each input: 

a 2 =a 2 a + a 2 + .... (17) 

Inputs to the fit include groups of priors associated with 
pieces of the fit function as well as statistical errors on 
the data points. Here we streamline the process by calcu- 
lating explicitly the differential of x 2 with respect to the 
inputs and so determining a 2 , a 2 etc. directly. The re- 



FIG. 7: Results for the mass of the D s meson tuned to the 
correct valence c and s mass on each ensemble from Table llVl 
as a function of the difference between the sea light quark 
mass and the physical value scaled by the physical strange 
quark mass (i.e. the parameter Sxi). The results are clearly 
separated by their lattice spacing value with very coarse at 
the top and ultrafine at the bottom. The lines show the result 
of the fit described by equation |15| taken at the value of the 
sea strange quark mass (Sx s ) and using the lattice spacing 
value corresponding to the ensemble with smallest Sxi in that 
group. The results on the coarse lattices at Sxi = 0.25 in- 
clude numbers at two different values of Sx a as well as at two 
different volumes. This gives an idea of the spread in results 
from these effects. The lowest line is the fit curve in Sxi at 
a — and Sx s — 0. The shaded band gives our final result 
adjusted for electromagnetic effects and with the full error as 
described in the text. The black burst gives the experimental 
result. 



suiting breakdown of errors given in Table [V] shows them 
to be dominated by statistical errors. 

Additional errors to be included in the error budget 
are errors that affect the final result in physical units 
but do not affect the fit above. The first of these is the 
overall error in the physical value of r\ of 0.7%. This 
affects the tuning of all the valence masses but, as de- 
scribed earlier, the effect on A is reduced by a factor of 
3 because of cancellation between scale shifts and tuning 
shifts. More precisely we find a 1.0 MeV error on 
from the r\ uncertainty. The effect on A of the uncer- 
tainty in the physical values of the rj c mass and the rj s 
mass used in tuning can be judged from Figures [3] and [4] 
The error on the r\ c mass has negligible effect, again be- 
cause most of the r\ c mass dependence cancels out in A. 
The uncertainty in the ?y s mass is not negligible, however, 
but gives an uncertainty in A, which we then transfer to 
jtod,, of 1.1 MeV. The error on the physical value of the 
r\ c reappears when we reconstruct m d 3 from A and m r)c . 
It therefore gives a 1.5 MeV error to mjj 3 coming from 
electromagnetic and annihilation effects in the r\ c meson 
mass. The error from electromagnetic effects on the D s 
mass itself is 0.7 MeV as described earlier. 
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TABLE V: Full error budget for mo s , fn s and f Vc given as 
a percentage of the final fitted value. Note that in the case of 
/,,,. the top six errors are those to be considered for a lattice 
QCD calculation that matches this one. As discussed in the 
text, the bottom three errors are included for completeness. 



B. 



h 



Error 




}d s 


fvc 


statistical/valence tuning 


0.094% 


0.57% 


0.45% 


n/a 


0.025% 


0.15% 


0.16% 


ri 


0.051% 


0.57% 


0.27% 


a 2 extrapoln 


0.044% 


0.40% 


0.24% 


m qySea extrapoln 


0.048% 


0.34% 


0.09% 


finite volume 


0% 


0.10% 


0% 


m vs 


0.056% 


0.13% 




em effects in D s 


0.036% 


0.10% 




em and annihln in m Vc 


0.076% 


0.00% 


0.05% 


em effects in rj c 






0.40% 


missing c in sea 


0.01% 


0% 


0.01% 


Total 


0.16% 


1.0% 


0.6% (top 6) 



The error from the finite volume of the lattices we esti- 
mate to be negligible from finite volume chiral perturba- 
tion theory. Our lattice results comparing two different 
volumes (sets 4 and 5) show no significant effect at the 
level of 0.4 

Our lattice calculation includes u, d and s quarks in the 
sea but no c quarks, although gluon field configurations 
are now being generated that do include them [39) . In 
the real world c quarks do appear in the sea and we can 
estimate the effect of these perturbatively because the c 
quark mass is relatively heavy, i.e. larger than typical 
momenta appearing inside the mesons we are discussing. 
The effect of a massive quark loop in the gluon propaga- 
tor which gives rise to the heavy quark potential is simply 
to add a correction to the potential which is proportional 
to a delta function at the origin [40 : 



V(r) = - 



C 



-c 



10mi 



5 3 (r) 



(18) 



Although this additional term is a spin-independent in- 
teraction its effects in charmonium can be judged by com- 
parison to that of the hyperfine potential. The hyperfme 
potential induces a mass splitting of « 120 MeV from a 
term which has the same 5 function form as above but a 
coefficient 280 (— 807r/(3a s ) ) times as large. Thus we 
expect the shift of the r\ c (and J/tp) masses caused by 
the presence of c quarks in the sea to be approximately 
0.4 MeV. The D s meson has much smaller momenta typ- 
ically inside it and so we expect a much smaller effect 
from c quarks in the sea on the D s meson mass. If we 
set that effect to zero, so that conservatively there is no 
cancellation of this effect in the quantity A, then we ob- 
tain an uncertainty in our final D s mass of 0.2 MeV, or 
0.01%. 

Our final result for mn s is then 1.9691(32) GeV to 
be compared to an experimental result of 1.9685(3) 
GeV HD. 



The decay constant of the D s meson is the main result 
from this paper. Having discussed in detail the tests that 
can be successfully done of the D s mass, we now discuss 
the analysis of the decay constant. 

Table [llT| gives the raw results for the decay constant on 
the 11 different ensembles we have studied. As for m,]j s 
it is important to be able to understand the dependence 
of fjj s on the valence c and s masses and to tune the 
result on each ensemble to the physical values for these 
masses. As described above, this corresponds to tuning 
them to physical values of the rj c and rj s meson masses. 
Figures [8] and [9] show the dependence of fo 3 on these 
meson masses on very coarse, coarse and fine lattices. 
Again we are using results somewhat above the physical 
values for the rj s and r\ c masses to extract the dependence 
which will then allow us to tune accurately our results 
that are much closer to the physical values. As expected, 
the dependence on = m s is linear and the slope 
does not change with lattice spacing. The value of the 
slope, 0.06GeV~ 1 can be compared to the change in fr> 
expected from q = s to q = I [3 . This corresponds 
to a somewhat larger slope of 0.09GeV _1 but is over a 
much larger range where nonlinear effects may appear. 
The slope of /n, against m ric falls from very coarse to 
fine lattices. This has interesting implications for the 
behaviour of the heavy-strange meson decay constant as 
a function of heavy quark mass. It is clear from the study 
of the T) s and D s mesons that the decay constant increases 
as the 'heavy' quark mass is increased from m s to m c . 
However, above m c the behaviour is less clear because 
lattice QCD calculations have so far not been accurate 
enough to distinguish clearly what is happening to within 
5-10% errors. There are known to be large corrections 
to the 1 / y/niQ behaviour expected from HQET because 
fn s and fs s are not very different [41]. This is consistent 
with a slope against heavy quark mass for fp s that tends 
to zero. It is clear that understanding this dependence 
also requires good control of discretisation errors. 

Again we use the dependence shown in these plots to 
make small tuning shifts to the values of fu s on each 
ensemble so that they correspond to the correct result 



for m„ c = 2.985 GeV and m Vs = 0.6858 GeV. Table IV 



gives the tuned values on each ensemble. Because the 
statistical errors are about twice as large for fo s as for 
A and the dependence on m, )c and m ri3 is smaller, the 
tuning shifts, and the errors from them, are very much 
less than the statistical errors on all ensembles. Even on 
set 2 the shift from mistuning is only la. The error in 
fu s from the uncertainty in r\ /a is only slightly reduced 
over its naive value from cancellations. It is also much 
smaller than the statistical error everywhere. It is given 
as the second error in Table IIVI 



Again it is clear from Table IV that the sea quark mass 
dependence of the results is smaller than our 1-2 MeV sta- 
tistical errors, but the lattice spacing dependence is not. 
We therefore fit the sea quark mass dependence with a 
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FIG. 8: Results for the decay constant of the the D a meson 
as a function of the square of the r] a meson mass, acting as a 
proxy for the strange quark mass. Results are for a range of 
different quark masses around the masses of the c and s quark 
masses on very coarse set 1 (triangles), coarse set 4 (crosses) 
and fine set 7 (squares). The lines are fits to the results for 
each ensemble allowing linear terms in m^ s and m J)t; . Here 
the lines join points for a fixed c quark mass. See Figure[9]for 
the equivalent as a function of m r>c . 



FIG. 9: Results for the decay constant of the the D s meson 
as a function of the r\ c meson mass, acting as a proxy for 
the charm quark mass. Results are for a range of different 
quark masses around the masses of the c and s quark masses 
on very coarse set 1 (triangles), coarse set 4 (crosses) and 
fine set 7 (squares). The lines are fits to the results for each 
ensemble allowing linear terms in m^ s and m, c . Here the 
lines join points for a fixed s quark mass. See Figure [8] for the 
equivalent as a function of m? a . 



relatively simple form that allows an error for what little 
dependence there is to be included in the final extrapo- 
lated value at the physical point. For the lattice spac- 
ing dependence we include relatively high order terms to 
make sure that a sufficiently large error is included in 
the final extrapolated value for this dependence. The fit 
form is the same as that used for A: 

4 

f Ds {a,Sx h 5x s ) = f Ds , P hys[l + ^2c j (m c a) 2: > (19) 

3=1 

+ 2bi5xi(l + c b {m c a) 2 ) 
+ 2b s Sx s (1 + c s (m c a) 2 ) 
+ Abu{8xif + 2bi s 5x l 5x s + b ss {Sx s ) 2 }. 

We take the same prior values and widths as before ex- 
cept that for fD s , P hys we take to be 0.25(10). 

The extrapolated result at the physical point, fD s ,phys 
is 0.2480(19) GeV with a x 2 /dof of 0.2 for 11 degrees 
of freedom. The fit is robust to changes in the fitting 
function: 

• changing the prior on all the Cj (including c\) to 
0.0(5) changes fD s ,phys by 0.8er and increases the 
error by 30%. 

• adding or subtracting two powers of a 2 into the sum 



on j in equation 19 does not change fo s .phys or its 
error. 

• adding an extra power of discretisation errors into 
both the linear and quadratic sea-quark mass de- 
pendent terms makes no difference. 



• missing out the sea quark mass dependence alto- 
gether changes fD 3 ,phys by 0.2ct but increases the 
X 2 value to 0.3. 

• Changing all the 5x values by 10% in either di- 
rection makes no appreciable difference, nor does 
changing them within their error bars on, for ex- 
ample, the ultrafine or fine lattices. 

• missing out the very coarse lattice results does not 
change fo s ,phys] missing out the very coarse and 
the coarse shifts fD 3 ,phys by 0.3cr (1 McV). 

• missing out the ultrafine result shifts fD s ,phys by 
0.4cr (1 MeV). 

Figure [10] shows the results plotted against the square 
of the lattice spacing. The line is the fit curve for the 
physical sea quark mass values (i.e. 6x\ = 5x s = 0). The 
shaded band is then the final physical result including 
the full error of 1.0% (2.5 MeV), to be discussed below 
and broken down into its component parts in Table |V| 

We construct the error budget as before, separating 
the error of 1.9 MeV resulting from the extrapolation 
to the physical point into its components of statistical 
error, r\/a error and errors from extrapolation in the 
lattice spacing and in the sea quark masses. Here the 
contributions from statistical errors and the different ex- 
trapolation errors are comparable. 

The error in the physical value of ri is 0.7%. This 
becomes a 0.6% error in fo s when the effects of r\ on 
shifting the value of m Vs are taken into account. The 
effect of the 0.6% uncertainty in the physical value of 
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FIG. 10: Results for the D 3 decay constant tuned to the 
correct c and s mass on each ensemble as a function of the 
square of the lattice spacing. The line shows the result of the 
fit at the physical value for the sea quark masses, as described 
in the text. The shaded band gives our final result with the 
full error bar as described in the text. 
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m rts can similarly be estimated from the dependence of 
fu 3 on the rj s mass at 0.1%. The uncertainty in fjj 3 from 
the uncertainty in the value of the rj c mass is negligible. 
The error from working on a finite spatial volume instead 
of infinite volume is estimated at 0.1% from comparing 
finite and infinite volume chiral perturbation theory. It 
is clear from our results (see Table III I that we see no 



significant volume dependence within our 0.5% statistical 
errors, which is in agreement with chiral perturbation 
theory, but that provides a stronger constraint. 

The size of electromagnetic effects inside the D s can 
be bounded by the size of these effects on the rj c . By al- 
lowing for an electromagnetic contribution to the heavy 
quark potential we estimate that f Vc could be increased 
by up to 0.4% by these effects. Since the D s has one 
quark of half the electromagnetic charge and is also much 
larger, so less sensitive to short-distance electromagnetic 
effects, we conservatively take an error of 0.1% from in- 
ternal electromagnetic effects [4"2"| . 

The error resulting from missing c quarks in the sea can 
also be bounded by the size of such effects on f, lc . In sec- 
tion [niX] we discussed a comparison between the hyper- 
fine potential in charmonium and that induced by adding 
c quarks in the sea. The hyperfine potential causes the 
difference between fj/$ and f Vc , which we will see in the 
next section is very small, 3%. The c-in-the-sea poten- 
tial is 280 times smaller and so will produce a completely 
negligible effect on f Vc and therefore also on fo s . 

Figure 11 shows the results for fo s as a function of 
the sea light quark mass, normalised to the strange mass 
as in equation |A3[ The lines show the fitted curves at 
the appropriate values of lattice spacing and sea strange 
quark mass, along with the final physical curve and final 
result with error band. No significant dependence on sea 
quark masses is seen. 

Our final result for f Ds is 0.2480(25) GeV, to be com- 
pared to the October 2010 average from the Heavy Flavor 
Averaging Group of 0.2573(53) GeV [25]. 
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FIG. 11: Results for the decay constant of the D s meson 
tuned to the correct valence c and s mass on each ensemble 
from Table IIVI as a function of the difference between the sea 
light quark mass and the physical value scaled by the physical 
strange quark mass (i.e. the parameter 5xi). The results 
are clearly separated by their lattice spacing value with very 
coarse at the top and ultrafine at the bottom. The lines show 
the result of the fit described by equation |19| taken at the 
value of the sea strange quark mass (5x s ) and using the lattice 
spacing value corresponding to the ensemble with smallest Sxi 
in that group. The results on the coarse lattices at Sxi = 0.25 
include numbers at two different values of Sx s as well as at 
two different volumes. This gives an idea of the spread in 
results from these effects. The lowest line is the fit curve in 
Sxi at a = and 8x s = 0. The shaded red band gives our 
final result with the full error as described in the text. 



Here we study the remaining independent quantity 
that can be extracted from the pseudoscalar correlators 
calculated here, the decay constant of the r\ c meson. Al- 
though this cannot be directly related to any process 
measurable in experiment, it can be compared between 
lattice QCD calculations using different formalisms for 
the c quarks. Since we have particularly accurate results 
here, we give a value for /„ that others can use to test 
their formalisms against. 

The raw results for f r/c on each ensemble are given in 
Table [TT] Since the r\ c contains only charm quarks we 
have only to plot f r)c against m 7)o to interpolate to the 
correct point on each ensemble. Because this is simpler 
than having to separate the dependence on two masses, 
as was done for the D s , we can plot the results from 



many more of the ensembles. Figure 12 shows the results. 



16 




3.1 3.2 3.3 

m„ (GeV) 



3.5 



FIG. 12: Results for the r] c decay constant as a function of 
the r)c mass for the different ensembles in Table [I] As in 
Fig. [I] very coarse ensembles are triangles; coarse, crosses; 
fine, squares; superfine, pluses. Errors shown are statistical 
only. Results for heavier sea u/d quark masses at each lattice 
spacing are in red, lighter ones are in blue. On the coarse 
lattices the very heavy sea masses of set 6 are in pink, the 
lighter strange sea mass of set 7 in grey and the large volume 
results on set 5 are in green, on top of the result from set 
4. The lines are fits to the results for one ensemble at each 
lattice spacing allowing linear and quadratic terms in m Vc . 



As expected, the dependence is linear (we allowed for 
quadratic terms in the fit, but these were small) but with 
a slope that depends on the lattice spacing. The figure 
also emphasises how little sea quark mass dependence 
there is, in line with the evidence from Figure [I] Some 
is visible above our very small statistical errors on the 
coarse and fine ensembles, however. 

Again we use the dependence shown in Figure [12] to 
make small tuning shifts to the values of f ric on each 
ensemble so that they correspond to the correct result for 



2.985 GeV. Table IV gives these tuned values. The 



statistical/tuning errors are small but the T\ja errors are 
even smaller because of cancellation when the retuning 
is done on changing the lattice spacing. Once again the 
lattice spacing dependence is the most striking feature of 
these results. 

We fit the tuned values to the same functional form as 
used for mo, (equation 15) and fo, (equation 19). We 



take the same prior values and widths for the parameters 
except that for the physical value of /„ o , fn ctP hys we take 
0.4(2) and for the coefficients, c,, for the discretisation 
errors we take 0.0(3), since v 2 for a c quark is expected 
to be somewhat higher than in a D s . 

The extrapolated value at the physical point, f rh , P hys, 
is 0.3947(20) GeV with a x 2 / dof of 0.3 for 11 degrees of 
freedom. Once again we tested how robust the fit was: 

• changing the prior on all the Ci (including c\) to 
0.0(8) changes fn c , P hys by 0.5er (1 MeV) and in- 
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FIG. 13: Results for the r\ c decay constant tuned to the correct 
c mass on each ensemble as a function of the square of the 
lattice spacing. The line shows the result of the fit described 
in the text. The shaded band gives our final result with the 
full error bar as described in the text. 
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FIG. 14: Results for the decay constant of the rj c meson tuned 
to the correct valence c mass on each ensemble from Table llVl 
as a function of the difference between the sea light quark 
mass and the physical value scaled by the physical strange 
quark mass (i.e. the parameter Sx{). The results are clearly 
separated by their lattice spacing value with very coarse at 
the top and ultrafine at the bottom. The lines show the re- 
sult of the fit described in the text, taken at the value of the 
sea strange quark mass (Sx 3 ) and using the lattice spacing 
value corresponding to the ensemble with smallest Sxi in that 
group. The results on the coarse lattices at Sxi = 0.25 in- 
clude numbers at two different values of Sx a as well as at two 
different volumes. This gives an idea of the spread in results 
from these effects. The lowest line is the fit curve in Sxi at 
a = and Sx s = 0. The shaded blue band gives our final 
result with the full error as described in the text. 



17 



creases the error by 40%. 

• adding two powers of a 2 into the sum on j in the 
fit equation does not change fn c ,phys or its error; 
subtracting two powers changes fr),.,phys by 0.5ct 
(lMeV) and reduces the error by 30%. 

• adding an extra power of discretisation errors into 
both the linear and quadratic sea-quark mass de- 
pendent terms makes no difference. 

• missing out the sea quark mass dependence alto- 
gether does not change f Vc , P hys but increases the 
X 2 value to 1. 

• Changing all the 5x values by 10% in either di- 
rection makes no appreciable difference, nor does 
changing them within their error bars on, for ex- 
ample, the ultrafine or fine lattices. 

• missing out the very coarse lattice results does not 
change f n ,.,phys appreciably; neither does missing 
out the very coarse and the coarse but the error 
increases by 50%. 

• missing out the ultrafine result shifts f Vo , P hys by 
1.4<t (2.5 MeV) and increases the error by 40%. 

The error budget is constructed as before, estimating 
the split in the error obtained from the fit into compo- 
nents from statistics, r\ja and extrapolations in a 2 and 
the sea quark masses. In addition the error from the un- 
certainty in the physical value of r\ becomes 0.3%, allow- 
ing for the cancellation that reduces the sensitivity below 
the naive 0.7%. The error from finite volume effects we 
take to be negligible based on the chiral perturbation 
theory studies of the much larger D s meson. 

f, la is not a quantity 



As we will discuss in section IV 



that can be compared directly to experiment. We include 
it here as a calibration point for lattice QCD studies of 
charm physics. As such, we do not have to include errors 
arising from effects outside a pure lattice QCD calcula- 
tion including u, d, and s sea quarks and taking the rj c 
mass to be 2.985 GeV. Thus in Table |V| only the top six 
errors in the final column should be included for such a 
calculation and the bottom three ignored. 

For completeness we discuss other sources of error that 
may need to be considered if lattice QCD calculations dif- 
fering in detail from ours are compared to it. The error 
that arises from the 3 MeV uncertainty in the physical 
value of the r\ c mass can be estimated from the slope of 

This gives an error of 0.05% 
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f Vc with m Vc in Fig 

with / r)e increasing with the value of . Internal elec- 
tromagnetic effects inside the rj c will also increase f Va . In 
section IIIB we estimated this effect at 0.4% (but lattice 



QCD calculations will not typically include electromag- 
netic effects) . The effect of including c quarks in the sea 
will also be to increase f Vc . In section |III B|we estimated 



Figure 13 shows against a 2 in fm 2 with the fit curve 
for the physical sea quark mass values. The shaded band 
is the final physical result including the full 0.6% error 
i.e. 0.3947(24) GeV. 



Figure 14 shows the results for f r)c as a function of 
the sea light quark mass, normalised to the strange mass 
The lines show the fitted curves at 



A3 



this as 0.01%, based on a comparison to fj/^ that will 
be described in section llV B 31 



as m equation 
the appropriate values of lattice spacing and sea strange 
quark mass, along with the final physical curve and final 
result with error band. No significant dependence on sea 
quark masses is seen. 



IV. DISCUSSION 

A summary of the results from this calculation is then: 

m Ds = 1.9691(32)GeV (20) 
f Ds = 0.2480(25)GeV 
f nc = 0.3947(24)GeV 

A. Comparison to our previous results 

Our new results improve on our 2007 results [5] in sev- 
eral ways, as described earlier. It is worth discussing 
the effect of these changes on the final numbers because, 
particularly in the case of fn e , the shift from 2007 is 
significant. 

Our 2007 result for m Ds was 1.962(6) GeV obtained 
from very coarse, coarse and fine ensembles. The lattice 
spacing was fixed using the quantity r\ as here, but set- 
ting the physical value of r\ to 0.321(5) fm. The error 
on tod, from this uncertainty in r\ was 0.2% i.e. 4 MeV. 
Since then we have improved significantly the calibration 
of the lattice spacing by improving the determination of 
the physical value of r\ to 0.3133(23) fm. This has used 
improved determinations of T\ /a on each ensemble by the 
MILC collaboration [57] . The change in the value of r\ 
represents 1.5cr and therefore we expect mu s to change 
by approximately 6 MeV. In fact the change has been 7 
MeV. 

Of course the change in n has not been the only 
change. The lattice spacing values on individual ensem- 
bles have moved relative each other with changes in T\ /a 
values. These have moved furthest on the very coarse 
set 1, changing by 1% or 2a, but with some changes of 
up to 0.5% (lfi) on the coarse ensembles. Values on the 
fine ensembles have not shifted signficantly. The relative 
shifts change the lattice spacing extrapolation slightly, as 
does our improved tuning of the charm quark mass (the 
strange mass was tuned within the chiral extrapolation 
previously using results for the K and 7r meson masses). 
We also have additional determinations of the sea quark 
mass dependence. These other effects largely cancel each 
other, however, in this case. Our new error budget shows 
an improved error coming from the determination of r\ 
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and this is the main effect behind the reduction of total 
error from 0.3% to 0.2%. 

Our f Dg result has changed by 2.7% (from 0.2415(32) 
GeV) which is a shift of 2a. From our error budget the 
change expected from the change in r% is 1.5%. Combined 
with changes in r±/a and improved tuning, however, re- 
sults on the fine and very coarse ensembles have changed 
by up to 2%. This has affected the continuum extrapo- 
lation. Sea quark mass effects, although not significant 
either now or before, have also changed in the same di- 
rection. This has meant that the 0.3% sea quark mass 
extrapolation error has added linearly to (some of) the 
0.5% continuum extrapolation and the roughly 2% shift, 
rather than in quadrature. 

The ratio JdJ Id is not very sensitive to r\ and so, 
although we have not yet performed an improved analysis 
of fo, we would not expect this ratio to change very 
much. If we take our previous result for /d s //d, but 
double the r\ uncertainty and add it linearly to the a 2 
and m Ui d extrapolation errors to allow for the behaviour 
seen in fo s we would obtain an error of 1.5% on the ratio, 
giving 1.164(18). Combined with our new result for fu g 
this gives a value for f D of 0.213(4) GeV, to be compared 
to a CLEO result of 0.206(9) GeV [13]. We emphasise 
that our new value for fz> does not result from a new 
analysis of fn itself but simply from the change in /rj a 
given here. 

Our 2007 results for and Jk change a little when 
the new value for is used. Using the fitting procedure 
described in the appendices of [34 (but not including 
the experimental results for and fx in the fit data), 
we find /„ = 132(2) MeV and f K = 159(2) MeV which 
agree within errors with our 2007 results [3] and with 
experiment [35] to within about 1.5<7. 



B. 



Comparison to other lattice results and to 
experiment 



mo. 



As discussed in section [ill A[ the accurate determina- 
tion of the mass of the D s meson is an important test of 
the calculation of fn s . Our result, 1.9691(32) GeV, is in 
good agreement with experiment, as shown in Figure [6] 
The experimental error is 0.3 MeV [35]. To improve the 
lattice QCD error of 3 MeV further would require im- 
proved statistical errors on the very fine lattices but also 
improved errors from elect romagnetic/?7 c annihilation ef- 
fects that are not currently included in lattice QCD cal- 
culations. It is impressive that lattice QCD calculations 
have reached the point where electromagnetic effects have 
to be considered in the match to experiment. 

Other lattice QCD formalisms for c quarks are not as 
highly improved as HISQ. They then have more difficulty 
in handling charmonium and so fix the c quark mass from 
the D s . However, we believe that it is still important to 
check the masses of other mesons containing c quarks as a 
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FIG. 15: Summary of full lattice QCD results for the dif- 
ference of binding energy between charm-strange and char- 
monium states. A uses pseudoscalar mesons D s and r/ c and 
compares the result from this paper to experiment, given by 
the appropriate red point and shaded band). A' uses a spin- 
average of the pseudoscalar and vector states and compares 
the result from the Fermilab Lattice/MILC collaborations to 
experiment. Our result corresponds to the complete error 
budget given in Table [V] and is corrected for missing electro- 
magnetic effects. The Fermilab Lattice/MILC result includes 
both errors given in [Hj but has not been corrected for miss- 
ing electromagnetic effects. 



test of systematic errors. The easiest quantity to compare 
is the one defined earlier as A = mo 3 — rn r)c /2, the dif- 
ference in binding energy between charmonium and D s . 
Our result for this is plotted in Figure [15] A variant of A 
was recently calculated using the Fermilab heavy quark 
formalism for c quarks, combining this with light asqtad 
quarks on the MILC very coarse, coarse and fine ensem- 
bles |43j . The c mass is fixed from the energy-momentum 
relation for the D s meson (because the energy at zero mo- 
mentum is not equal to the mass), which leads to size- 
able statistical errors in the tuning process, growing with 
heavy quark mass jH]. Typically the 'kinetic mass' for 
the D s is obtained to 2%. The Fermilab lattice/MILC 
collaborations quote a result for A' = m(D s ) — m(lS)/2 
of 0.529±7+ 1 Q 2 with a partial error budget g3]. Here TT S 
indicates the spin-average mass of the D s and the D* 
and m(lS) is the spin average of the masses of the J/tp 
and the r\ c . The spin average is used to reduce their dis- 
cretisation error from spin-dependent terms, but the D* 
does have a strong decay mode, albeit Zweig-suppressed, 
that will lead to an additional systematic error in the 
lattice QCD calculation. The first error given above is 
from statistics and extrapolation uncertainties and the 
second from the physical value of r\ which they take 
as 0.318 



+0.000 
-0.007 



fm. The Fermilab Lattice/MILC result 
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FIG. 16: Comparison of the result from this paper for the 
D s decay constant with those from other lattice QCD calcu- 
lations that include the effect of sea quarks. The Fermilab 
Lattice/MILC result is a preliminary one but also includes 
the effects of u, d and s quarks in the sea. The ETMC re- 
sult includes only u and d quarks in the sea. We show also a 
recent average of experimental results from the Heavy Flavor 
Averaging Group [25] and two separate averages over the 
and rv channels. Experimental results for fn s convert the 
lcptonic decay rate to a decay constant using equation [21] and 
using an input value for V C3 (see text). 



agrees with experiment and is plotted in Figure 15 for 
comparison to our result for A. More detailed com- 
parison between the results needs improved accuracy for 
those from the Fermilab formalism. 



Figure [16] compares the result for the D s decay con- 
stant from this paper to other lattice QCD calculations 
that include the effect of sea quarks. The Fermilab Lat- 
tice/MILC result of 260(10) MeV is a preliminary one 
from a conference presentation |45) . updated from their 
original 2005 calculation [7] in a number of ways but in- 
cluding an update of the physical value of the parameter 
T\ used to set the lattice spacing as we have done here. 
Their calculation uses MILC gluon field configurations 
as we do, but at the three coarsest lattice spacing values 
that we have used. The Fermilab formalism for c quarks 
is combined with the asqtad formalism for the s quarks. 
As explained above the c quark mass is tuned from the 
dispersion relation for D s mesons. In the Fermilab for- 
malism there is no PCAC relation and so the temporal 
axial current operator that annihilates the D s in its lep- 
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FIG. 17: Comparison of our new result for the D s decay con- 
stant with recent experimental results from CLEO [21-23] and 
BaBar [26]. These are derived from leptonic decay modes of 
the D s in various channels, and using equation [2l] with an 
input value for V cs (see text). The CLEO numbers are taken 
from the compilation in [53], using consistent values for V C3 , 
uld b and td„ and so differ slightly from the historical numbers 
in Figure [18] We also include the HFAG 2010 world average 
for experiment |25| . 



tonic decay (equation 11 ) must be renormalised to match 
the continuum current operator that couples to the W. 
This is done by a perturbative calculation to 0(a s ) af- 
ter taking a ratio to vector current operators. The sys- 
tematic uncertainty from this approach is in principle 
0(a 2 s ) (« 5%), but it is argued in [45l [46] that a signifi- 
cantly smaller (1.4% + 0.3%) error be used which is the 
square of the one-loop contribution. It would be useful 
to test this on a calculation such as fx where the result 
is well-known [47]. With relativistic formalisms such as 
the HISQ formalism used here and the twisted mass for- 
malism to be discussed below, the existence of the PCAC 
relation means that the issue of renormalisation does not 
arise. Also in both cases, fx can be calculated as well as 
fo as a test of the error analysis. 



Figure 16 also includes the result 244(8) MeV from 
the European Twisted Mass Collaboration [?5] using the 
twisted mass formalism for all of the quarks. This for- 
malism is an improved version of the Wilson formalism 
with discretisation errors starting at 0(a 2 ), somewhat 
worse than the 0(a s a 2 ) for HISQ, but also having a par- 
tially conserved axial current so no renormalisation is- 
sues. ETMC include only the effect of u and d quarks 
in the sea, however, and it is not clear what systematic 
error to take for missing s quarks that are there in the 
real world. We cannot use perturbative arguments, as 
we have done here to account for the missing c quarks in 
the sea. ETMC are now improving their calculations to 
include both s and c sea quarks [49J. 
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The experimental results shown on figure 16 are the 
October 2010 averages from the Heavy Flavor Averaging 
Group [25], using recent CLEO [HH23], BaBar [H d] 
and Belle [TS] results from measurement of the D s — » /.iv 
and D s — > tv decay rates. To determine fo 3 from experi- 
ment the measured leptonic branching fraction, corrected 
for electromagnetic radiation [32], is used in: 



Jd s = 



1 



G F \V cs \mi(l - mf/m 2 D ) 



I8itB{D s -> Iv) 
m D T Ds 



(21) 



A value for V cs must be assumed. In the past V cs = V uc i 
has often been taken (see, for example, (21]'). assuming 
2x2 CKM unitarity. HFAG take the 2010 Particle Data 
Tables result for V cs (0.97345(16)) from a full CKM ma- 
trix unitarity fit [251 135] . These two alternatives for V cs 
differ at the level of 0.1% which is irrelevant here. 

It is clear from Figure [16] that there is no longer any 
significant 1 fn B puzzle' [51] since the discrepancy between 
our lattice QCD result and the world average of exper- 
iment (257.3(5.3) MeV) is 1.6a. The average of exper- 
imental results in the tv channel (252.4(6.9) MeV) and 
our value agree very well. This is emphasised further in 
Figure [17] where the most accurate recent experimental 
results are individually compared to our value for fo s , 
and all except one disagree by less than la. 

Things have now changed quite significantly since the 
summer of 2008 when the most accurate experimental re- 
sult for fjj t was 267.9(9.1) MeV [H] and the most accu- 
rate lattice QCD result was 241(3) MeV [3] E], differing 
by almost 3 a. The experimental average moved down 
5% (1.5a) in early 2010 but has since moved up 1% to the 
new world average value and the lattice result has moved 
up 3% (2.3c). The discrepancy between experiment and 
lattice QCD is now only 4% (1.6 a) and the experimental 
error is now reduced to only twice that of the lattice QCD 
error. This marks significant effort both experimentally 
and theoretically on this quantity to understand and pin 
down the original discrepancy. Figure [18] shows the his- 
tory of fu s from experiment and lattice QCD since the 
first full lattice QCD calculation of 2005. 



As stated earlier, there is no direct comparison pos- 
sible between lattice results for f nc and experiment be- 
cause the r\ c does not annihilate to a IF boson or other 
particle that would couple directly to the temporal axial 
current. The high accuracy of our results is therefore use- 
ful only to provide a comparison point for other lattice 
QCD calculations. No result of comparable accuracy is 
available from any other charm quark formalism as yet 
and including the effect of sea quarks. ETMC [52] quote 
a preliminary result of 379(29) MeV for f Vc including u 
and d quarks in the sea only and tuning the c mass from 
the mass of the J/ip (i.e. this analysis is not directly 
linked to their D s analysis, as ours is). Future lattice 
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FIG. 18: Values for fr> s from experiment and from lattice 
QCD since 2005, excluding results from conference proceed- 
ings. Later results from a given collaboration and process 
supersede the earlier ones. Experimental results are divided 
into those from the \w channel 41 12, 15, 2 1 [26] (in red) and 
those from the TV channel [TT] Q31 |2"TH2"31 |2"o] (in several r de- 
cay modes, in pink). The HFAG October 2010 world average 
for experiment [25] is included as a light orange band. Note 
that the leftmost red point (from BaBar [T2]) appears with 
dashed error bars - the lower value with solid error bars is the 
result adjusted by HFAG [50], although this number is not 
now included in the HFAG average. Lattice QCD points are 
in dark blue for full QCD [3] [7] - the rightmost point is from 
this paper. The light blue point is from ETMC |48| including 
only u and d in the sea. 



charmonium calculations using different formalisms (for 
example [53] or [54] ) can use our result as a benchmark 
point to check renormalisation or discretisation effects 
because f Vc is a very simple quantity to calculate. 

Although direct comparisons with experiment do not 
exist, various comparisons that rely on approximation 
schemes, principally potential models, can be made. In 
a potential model the decay constant of an S-wave state 
is related to the wave-function at the origin, ^(0), by 
"0(0) = fy/M/12, where M is the meson mass and / 
its decay constant. This relationship is only correct up 
to relativistic and radiative corrections, which for the rj c 
could be sizeable (at the level of 30%) . Using this same 
potential model approach the leading term in the decay 
width for r\ c 77 can be written as [55] : 



r(?7c -> 77) 



12ne 4 c a 2 \iP{0)\' 



(22) 



Here the c quark has electromagnetic charge e c (in units 
of e) , mass m c and a is the electromagnetic coupling con- 
stant. This formula has radiative and relativistic correc- 
tions at the next order. The decay width is only poorly 
known for the rj c with the PDG estimate given as 7.2(2.1) 
keV [35] • Substituting the decay constant into the for- 
mula and taking m c = M I)c /2, justifiable at this order, 
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gives f Vc = 0.4(1) GeV, where only the large error from 
experiment is shown. Alternatively one may extract f ria 
from B decays to r\ c K using the factorization approxi- 
mation. CLEO obtain f Vc = 0.335(75) GeV |£B]. 

A more useful experimental result to compare to our 
decay constant is probably to the decay constant of the 
J/ip. Because the J/ip can annihilate to a photon (seen as 
two leptons in the final state) through the vector current 
there is an exact relationship between the decay width 
and the decay constant of the vector particle defined in 
an analogous way to that for the pseudoscalar meson by 



^2 < 0\^\Vi > /3 = f v m v . 



(23) 



This decay constant can also be calculated in lattice 
QCD [51] . Work is in progress and results will be given 
elsewhere. The relationship between decay width and 
decay constant for the process Vh — > e + e~ is then 



T(V h 



4tv 



~^~ a QED e Q 



m v 



(24) 



The experimental results [35] give fj/^ — 407(5) MeV 
using l/aQED(m c ) = 134 [ST]. Thus 1% accurate re- 
sults for this decay constant are available from experi- 
ment, and can be used to test lattice QCD. In a poten- 
tial model vector and pseudoscalar values of -0(0) should 
differ only by relativistic corrections, since this is a spin- 
dependent effect which appears first at sub-leading or- 
der in the velocity-squared of the heavy quark. Thus 
we would expect our results for the pseudoscalar decay 
constant to be fairly close to those for the vector. It 
is hard to make this statement quantitative however be- 
cause, even if the difference in values of vector and 
pseudoscalar were accurately pinned down, the relation- 
ship of ip(0) to the decay constant could have sizeable 
radiative and relativistic corrections. 

Our result for f Vc , 0.3947(24) GeV, is in fact very close 
to the experimental result for fj/^p, only differing by 3% 
(2er). This is somewhat surprising, given naive potential 
model arguments. Accurate lattice QCD studies in bot- 
tomonium will show whether this is a coincidence at the 
charm mass or a more general feature. 



V. CONCLUSIONS 

In this paper we have updated our 2007 result for the 
mass and decay constant of the D s meson [3] to incorpo- 
rate a new more accurate calibration of the energy scale 
in lattice QCD. We have also included results at two finer 
values of the lattice spacing so we now cover a range of 
lattice spacing values from 0.15 fm down to 0.044 fm for 
improved determination of the continuum limit. Our re- 
sults for m jj s and fo s increase as a result of this calibra- 
tion. mo 3 is in excellent agreement with experiment with 
a reduced (3 MeV) error to give 1.9691(32) GeV. Our re- 
sult for fo s has increased significantly to 0.2480(25) GeV. 
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FIG. 19: Bounds in the tan/3/charged Higgs mass plane 
for a 2-Higgs doublet model of Type II given by our lattice 
QCD (i.e. Standard Model) result for fn s and two different 
experimental scenarios. The green area is excluded at 3a by 
existing experimental results [55] and the light blue area will 
be excluded by BESIII results 58 if the central experimental 
value does not change. The light grey band shows the direct 
limit from LEP searches 1601. 



This, along with recent movement of the experimental 
results, means that the '/d, puzzle' is essentially solved: 
there is no longer significant disagreement between the- 
ory and experiment for this quantity. The experimental 
error is double the theoretical error, however, and im- 
proved experimental results from BESIII aim to obtain a 
1% on fo s [58 . The lattice QCD error could be further- 
reduced by improved statistical accuracy on the very fine 
lattices. 

Instead of assuming a value for V cs to obtain an exper- 
imental result for fjj s to compare to lattice QCD we can 
combine our result for fjj s with the experimental leptonic 
branching fraction to give a direct determination of V cs . 
To do this we take the HFAG determination [25] of the 
world average leptonic branching fractions for the D s to 
[iv and tv of 0.590(33)% and 5.29(28)% respectively, our 
result for fo and 



Vr, = 



1 



^ttB{D s -> Iv) 



G F f Ds mi(l - mf/m 2 D J ]] 

This gives results for V cs of : 

V cs = 1.033(31), D 
= 0.990(28) 



m D t d 



De — s- tv 



(25) 



(26) 



where the error is dominated by the experimental branch- 
ing fraction. We can combine the results, allowing for 
correlated errors in fp s and td s , to obtain 



V cs = 1.010(22) 



(27) 
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This central value is in a disallowed region above 1 so we 
also provide an alternative result that takes this into ac- 
count. We divide the error above into its statistical and 
systematic contributions as 1.010(20)(11) and then rein- 
terpret the statistical probability distribution as a Gaus- 
sian cut off at 1. We then take the central value as the 
median of this new distribution and the error bars as en- 
compassing ± one third of the area about the median. 
This procedure gives the following result: 



0.990 



-0.007 
-0.012 



±0.011. 



(28) 



Both these values for V cs are compatible with CKM re- 
sults (or V u d) at better than the 2a level. An independent 
direct determination of V cs is possible from D — > Klv 
semileptonic decay for which it is also possible to obtain 
very accurate results with the HISQ action [59] . 

A useful bound can be obtained on the mass of a 
charged Higgs from comparing the experimental deter- 
mination of the D s leptonic branching fraction to the 
expected result using fo s from lattice QCD (i.e Stan- 
dard Model), see, for example, [2D]. In a 2-Higgs doublet 
model (Type II) the D s can also annihilate to a charged 
Higgs which interferes destructively with the W annihi- 
lation. This changes the leptonic branching fraction by 
a simple factor r, where 



<r = 1 



1 



1 + m s /m c \m,H 



1 



tan 2 /3 (29) 



and tan j3 is the ratio of vacuum expectation values of 
the two scalar doublets, r < 1 for large tan/3 but this 
would be seen from an experimental determination of 
fr> s (using V cs from CKM unitarity) being smaller than 
the lattice QCD result. Thus we can derive a bound in 
the tan/3/m#± plane from the fact that this is not the 
case. Here we update what was done in [20] to include 
our new lattice QCD result given here and the current 
world average fo s from experiment |25j . These combine 
to give a central value and error for -J¥= 1.038(23), i.e. 

using our recent 



2') 



\Jr > 0.968 at the 3<r level. Equation ! 
accurate determination of m c /m s from lattice QCD 
then excludes low values of m H ± as indicated in Fig- 
ure 19 The bound is not as strong as in (2UJ because 
of the upward shift of our lattice QCD result. However 
the fact that our lattice result, and now the experimental 
average, are so accurate still means that a bound exists. 
New results from BES [58] with improved experimental 
errors would produce a much stronger bound, if the ex- 
perimental central value does not change but the error on 
fu B is reduced to 1%. This is also indicated in Figure 19 
The exclusion limits should be compared to that from 
direct searches at LEP (m H ± > 78.6 GeV at 95% C.L.) 
from [5DJ and the estimates of discovery potential and 
exclusion reach of ATLAS at LHC [51]. Ref [52] obtains 
a bound of m H ± > 316 GeV from combining results from 
several processes including D / D s leptonic decay. 

We have also updated results for f„, f K and fr> based 
on the change in the calibration of the lattice spacing 



used here for fu s but, however, with no new calculations 
in these cases. We find results consistent with experi- 
ment. Finally we have given a new very accurate result 
for / r)c which will be useful as a calibration point for fu- 
ture lattice QCD calculations in charm physics. 

Acknowledgements We are grateful to the MILC 
collaboration for the use of their configurations. Com- 
puting was done at the Ohio Supercomputer Center 
and the Argonne Leadership Computing Facility at Ar- 
gonne National Laboratory, supported by the Office of 
Science of the U.S. Department of Energy under Con- 
tract DOE- AC02-06CH1 1357. We acknowledge the use 
of Chroma 63 for part of our analysis. This work was 
supported by the Scottish Universities Physics Alliance, 
STFC, MICINN, NSF and DOE. 



Appendix A: Sea quark masses 

The staggered quarks in the sea are asqtad improved 
staggered quarks whereas the valence quarks are HISQ 
quarks, i.e. they use different discretisations of the quark 
piece of the QCD Lagrangian. The s quark mass in 
the two formalisms will then not be the same, but there 
should be a fixed ratio between the two which is in princi- 
ple calculable in perturbation theory up to discretisation 
effects. This reflects the fact that the difference between 
the two Lagrangians is a difference of regularisation and 
therefore an ultraviolet effect. Calculations in 0(a s ) per- 
turbation theory of the mass renormalisation in the two 
formalisms shows that the 0(a s ) term in the relative nor- 
malisation is very small [29, 53] . We therefore expect [55] 



am 



hisq 



am asq 



1 -0.004a s (a) + Ca 2 s (a) 



(Al) 



up to discretisation and sea quark mass effects. Here 
am hisq an( j am asq are foe lattice valence quark masses 
for the HISQ and asqtad actions respectively that give 
the same meson mass for a particular meson on a given 
ensemble. Note that am asq is defined in the conventional 
way i.e. without the Ug factor present in Table [T] 

Given the HISQ to asqtad mass ratio we can determine 
the tuning of the sea quark masses from our tuning of 
the valence HISQ masses. There is very little sea quark 
mass dependence in the quantities that we study here, 
so that we do not need to know this ratio accurately. 
In principle it could be done very accurately, because as 
we have seen the meson masses can be determined very 
accurately. In the absence of this information for asq- 
tad quarks, however, we take the suggested tuned asqtad 
strange quark masses from the MILC collaboration 27J, 
correcting for the uq factor (taken from the lightest sea 
quark mass ensemble at each lattice spacing and given in 
Table VI ), and compare them to our tuned HISQ strange 
quark masses [34] . Figure 20 shows results on very coarse, 



coarse, fine and superfine lattices. The errors on each 
point are substantial, ~ 3%, because we have included 
the tuning error from each action added in quadrature, 
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TABLE VI: The sea asqtad masses given in table [T] have a 
factor of Mo equal to the fourth root of the average plaquette 
included in them. We remove this factor in our comparison 
of quark masses between HISQ and asqtad and so give values 
here in column 2 from [BB] , for all ensembles except set 7 where 
the result is simply estimated from that of the other coarse 
lattices. Column 3 gives values for the physical asqtad strange 
quark at each lattice spacing mass quoted by MILC [27] and 
including the no factor. The result on set 11 is obtained from 
the tuned HISQ strange mass and the ratio described in the 
text. Columns 4 and 5 then give values for Sxi and Sx s as 
defined in equation | A3| and used in our extrapolations to the 
physical point. Errors come from the errors in uoam 3tP h y s and 
are correlated between ensembles at a given lattice spacing 
and between Sxi and 8x s 



Set 


MO 


asq 


Sxi 


8x s 


1 


0.8604 


0.0439(18) 


0.184(10) 


0.10(5) 


2 


0.8610 


0.0439(18) 


0.405(19) 


0.10(5) 


3 


0.8678 


0.0350(7) 


0.106(3) 


0.429(29) 


4 


0.8677 


0.0350(7) 


0.249(6) 


0.429(29) 


5 


0.8677 


0.0350(7) 


0.249(6) 


0.429(29) 


6 


0.8688 


0.0350(7) 


0.535(12) 


0.429(29) 


7 


0.868 


0.0350(7) 


0.249(6) 


-0.143(18) 


8 


0.8782 


0.0261(5) 


0.201(5) 


0.188(23) 


9 


0.8788 


0.0261(5) 


0.439(10) 


0.188(23) 


10 


0.8879 


0.0186(4) 


0.157(5) 


-0.03(2) 


11 


0.8951 


0.0135(5) 


0.170(8) 


0.04(4) 



since the tunings were done in a different way. The re- 
sults can easily be fit to the form: 



am 



hisq 



am asq 



= 1 - 0.004ev s (a) + Ca 2 s (a) + Da 2 + Ea 4 (A2) 



The 



adding discretisation errors to that in equation |A1| 
fit gives a coefficient C«2. 

The fitted curve enables us to determine that the sea 
strange quark mass on the ultrafine lattices should be 
0.0135(5) (with uq factor included) i.e it is reasonably 
well-tuned. The error is substantial, but the sea quark 
masses have very little impact on the accuracy of results 
given here. 

gives values for the uq parameter (—(plaq) 1 ^ 4 ) 



Table 



VI 



and the physical asqtad strange quark masses given by 
the MILC collaboration [57] for all lattice spacing values 
except the ultrafine. We use our result for ultrafine as 
discussed above. From the physical strange quark mass 
we determine the physical light quark mass using the 
MILC result: m s /m ; = 27.2(3). The table then gives 
values for 5xi and Sx s , where 



0X q = 



^q,sea,phys 



^s,sea,phys 



(A3) 



used in our extrapolation to physical quark masses in 
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FIG. 20: Results for the ratio of the physical masses of the 
strange quark using the HISQ formalism [31] and the asqtad 
formalism |27| . The mass for the asqtad formalism in this ra- 
tio has the factor of uo removed. The points are lattice results 
combining our work and that of the MILC collaboration. The 
shaded band represents a fit of the form given in equation|A2| 
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